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We derive two new equations of quantum gravity and combine them with reinter- 
pretations of previously proposed concepts of dark energy, black holes, inflation, the 
arrow of time and the characteristic energy at which rest-mass first manifests itself 
into a theory which may be a first step toward a comprehensive description of all 
these phenomena. The resulting theory also predicts new tests which can be experi- 
mentally checked within just a few years. The two new equations are : A) a creation 
equation to give stimulated emission for any surface filled with gravitons, pulling 
energy from a background, and B) the association of an outgoing soliton wave of 
gravitons, a "shell front" with a large Lorentz factor derived from the uncertainties 
in both space and time. We model the background as a strong gravity brane, a 
Planck length apart from our brane in a fifth dimension. These new equations are 
combined with the common notions of an all-pervasive background of gravitons at 
the Planck limit, the "Planck sea"; the identification of the thermodynamic limit 
with the emission of gravitons in a "shell front" , i.e. what is usually called the en- 
tropy of black holes is identified with the outgoing gravitons; the concept of black 
holes as a membrane full of gravitons at a large Lorentz factor, the "Planck shell" ; 
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the emission of gravitons created in a "horizon shell" during inflation, which lose 
energy adiabatically with the Planck time. These equations result in stimulated 
emission of gravitons by interaction with the background, the "Planck sea" , to de- 
scribe dark energy, black holes, and the inflationary period of the universe. These 
proposals lead to gravitational waves constituting dark energy. These waves should 
be detectable within a few years with pulsar timing arrays. The extremely high, 
but finite Lorentz factor for signal propagation may be expected to have further 
consequences in particle interactions. 

PACS numbers: 98.80.Bp, O4.60.Bc, 95.36,+x 



I. INTRODUCTION 

The goal of this paper is to provide the first step towards a description of a compre- 
hensive model of the birth and evolution of the universe. This model describes how the 
interaction with an all-pervasive background field drives inflation, black holes, dark en- 
ergy, and the arrow of time. The background field is assumed to be the always the same 
one, sending energy into our world by stimulated emission. For a review of the standard 
quantum mechanical picture of the evolution of the universe, see pp. 

One of the new equations is the creation equation derived in section III. The other 
is the Lorentz factor expression for the motion of the shell front of emitted gravitons 
derived in section IV. The form of the creation equation is similar to the one obtained 
from the non-relativistic form of the Boltzmann equation (the Kompaneets equation), and 
the square of the phase space density occurs as in both the relativistic and non-relativistic 
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forms of the Boltzmann equation. The argument for the Lorentz factor is similar to the 
connection between the scale of the universe with the Planck scales as done in [2 |3], as 
well as the equation describing the time delay of arrival between a photon and an ultra 
high energy particle emitted simultaneously from across the universe. Our concept is quite 
specific, and leads to predictions, which can be either verified or disproved within a few 
years, probably within existing experimental plans. We generally use a micro-canonical 
ensemble approach, a statistics method (e.g. [HE]), instead of thermodynamics. It is well 
known that gravitational systems are never in full thermodynamic equilibrium, and so for 
gravitational systems the expression giving no change of entropy 

dS = (1) 

never holds. The microcanonical ensemble is thus the appropriate ensemble to use for the 
description of gravitational systems. 



The creation equation (37) is a combination of the "holographic principle" [HIE] and 
the relativistic form of the Kompaneets equation, [5HTU]. In our model a quasi static back- 
ground field exists which can emit gravitons via stimulated emission into our universe. 
The connection between the background field and our world is in "shells" whose thick- 
nesses equal the Planck scale ~ 10~ 33 cm in the comoving frame. The thickness of a shell 
is determined by the size of the fifth dimension. This background field and the emitted 
gravitons are taken to be real. The background field will be referred to as the "Planck 
sea", and it is described in the observer frame. Its energy spectrum is a Planck spectrum 
at maximal energy density, peaking at the Planck energy as seen by an observer in our 
world. This corresponds rather closely to the notion of an "all-pervading background 
field", used quite often in the literature ( e.g. [H]). The distinction between our model 
and what is described in [TTJ is that we assume that the background is important at all 
times. We will show below how the source of energy in our universe can be incorporated 
into the energy-momentum tensor and the Einstein equations, and the form is well known. 



4 



The approach described here may also explain the inflationary period, with new gravi- 
tons produced in a horizon shell, provided that the newly created gravitons lose their 
energy adiabatically with a time scale proportional to the Planck time, independent of 
redshift within the inflationary period; this is just a specific application of standard mod- 
els. In this concept the universe starts with one Planck energy graviton in one Planck 
volume; this is conceptually similar to some ideas of Lemaitre |12J. In our model the 
inflationary period ends when the number density distribution of the created gravitons 
begins to deviate significantly from the Bose-Einstein distribution, which coincides with 
the first appearance of the Higgs boson. 

The gravitational waves constituting dark energy in our concept may be measurable 
by pulsar timing arrays. The continuous creation of gravitons due to the interaction of 
Planck shells with the background implies an energy density - pressure relation 

Pde = Pde c 2 . (2) 

The gravitational waves produced in the formation of the first generation of stellar black 
holes should be measurable by gravitational wave interferometers. The soliton-like passage 
of a single graviton shell front may be measurable by Lunar laser ranging, as well as by 
ultra-precise timing measurements. 

This paper is organized as follows: In section II we give a detailed description of the 
background Planck sea. In section III we derive the creation equation. In section IV we 
discuss dark energy, in section V black holes, in section VI inflation, in section VII predic- 
tions for possible observations, and in section VIII we give some further consequences such 
as an emergent arrow of time, and the first appearance of rest mass, as well as summarize. 
For numerical illustration we use black holes at 3 10 6 M and a formation redshift of 50. 
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II. THE PLANCK SEA 

We interpret the usual notion of a quasi static background field as one filled with 
gravitons. Models with a static background have been treated in the literature, see for 
example [13J. Strictly speaking the background is evolving on a time scale which is slower 
by the energy density ratio, = 10 115 - 3 erg/cm 3 in the background (the Planck limit) to 
€de — 10~ 8 erg/cm 3 in our universe (the dark energy) as we will show below in section 
IV. C and VI. A possible model which has the features we require is two four-branes 
embedded in a five-dimensional bulk and separated by the Planck length; such a model 
is inspired by but different from the Randall-Sundrum models [HI [15]. A metric which 
describes the transfer of energy from a four-brane where gravity is strong to the brane on 
which we live is given by 

ds 2 = _ e (u/irt/^ c 2 dt 2 + e (l- b (H)" )2t/QrHe _ (f) fc (l ^)^^ +e(frt/ ^ M2 (3) 

where i — 1,2, 3, th is the Hubble time, I is the Planck length, u is the coordinate in the 
fifth dimension, and the remaining, non-coordinate quantities are arbitrary parameters. 
This metric is not an exact solution of Einstein's equations since the five-dimensional 
covariant divergence of the energy-momentum tensor does not vanish everywhere. How- 
ever, the five- dimensional divergence vanishes on the weak-gravity brane (u = 0), and 
it is approximately zero on the strong-gravity brane (u = I) for large (3/th and ip/r^. 
This metric describes a weak-gravity brane which is expanding with time and a strong- 
gravity brane which is contracting with time, albeit very slowly for the latter brane. The 
cosmological constant measured on the weak-gravity brane is 

^-weak— gravity 3 - — (4) 

The evolution of the strong-gravity brane is determined by the relative sizes of the 
parameters a , r# and (3. The metric in Eq.([3]) bears a resemblance to the first Randall - 
Sundrum model except that the metric tensor elements in Eq.([3l are time-dependent, and 
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the tensor element #44 depends exponentially on the 5th coordinate, u p3lfT5]. This metric 
is clearly an approximate description of the current geometry of the universe and does not 
describe the geometry of the universe for all epochs. The vanishing of terms linear and 
quadratic in time appearing in the Einstein tensor, Goo, evaluated on the strong-gravity 
brane requires that 

kar H 2kar H (3 
2b n A [3 no + par H 

and, using the first expression for above, 

k = , or k = — — . (6) 

A(3nb + par H w 

The scale for changes in time is given by the ratio 

^change Tff (7) 

The scale of changes in time is thus, including factors of (4-7r) 3 , 

Tchange ~ 10 12 ° T H . (8) 

Factors of this magnitude recur often in cosmology, as already noted by Dirac [TSj, who 
showed that the mass of the universe is about 10 80 the mass of a proton, while other 
factors such as the ratio of the electric to the gravitational force are on the order of 
10 40 . The factor here is approximately 10 120 [El HE], and the long-sought explanation for 
this large factor emerges naturally from our model. We call this field the "Planck sea", 
and assume that the gravitons in this sea have a Planck spectrum at maximal energy 
density peaked at the Planck energy. The Planck sea is in equilibrium with respect to 
gravitational collapse, because the free-fall time scale and pressure wave time scale are 
equal at the Planck length. Larger scales are not causally connected. As noted in the 
Introduction, the background field and the emitted gravitons are not virtual but are taken 
to be real. The background field spectrum is modified at the (lower) energy where the 



7 



coupling strength of the four natural forces diverge, about a factor of order 10 18 from 
the maximum in the very early universe, as suggested by observations of the end of the 
inflationary period. This field provides the energy for the gravitons created by stimulated 
emission in our universe, as we discuss below. The total of the background field energy 
and the energy in our observed universe is conserved. 
We assume a specific field with a spectrum of 

vr 2 e e - 1 

in the observer frame where 



r2 .e ( 9 ) 



e = 2 10 ) 

mpic 2 

giving maximal density and maximal temperature, where ooqw is the frequency of the 
gravitational waves, and the other symbols have their normal meaning (see section IX.F.l, 
List of definitions). This spectrum holds down to those energies where the number density 
of resonant gravitons decreases even more steeply due to the separation of gravity from 
the other natural three forces, as one can show by solving the Boltzmann or Kompaneets 
equation on an FRW universe. 

In four space-time dimensions the relativistic Boltzmann equation on the weak-gravity 
brane (for a detailed discussion see Appendix A) can be written in terms of a parameter 
x which for convenience is defined as 

kb Jo 

where vq is the frequency of a wave at emission, and To is a temperature which is char- 
acteristic of the background source of energy. In terms of x the expression for ftf(v ,t) 
is 

>™ = £k%m"™- (12) 



s 



where T g o is the initial graviton temperature, and kq is the initial phase space integral of 
the matrix element squared for quadrupole emission of a graviton wave from the back- 
ground. 

This expression can be put into a more compact form by defining the dimensionless 
parameter, y, 

K B 1 1 



or 



T T 

This gives a spectrum of, (compare with Eq. [9 



y=l t ^r^dt> (14) 



M(x,y) = ±^- (15) 
n z e A — 1 

with 

x = yp+M (i6) 

where y is a time- dependent term, not significant until the square of the dimensionless 
frequency x matches y, so in the late universe the spectrum is steepened. As stated in 
the Introduction, this steepening coincides with the appearance of the Higgs boson. For 
a detailed derivation of Af(x, y) see Appendix A. 

Invoking a resonance at energies below this implies a reduced density, and so will 
violate the surface density condition in the early universe, at the end of inflation. Any 
attempt to resonate with energies beyond the Planck energy also fails. The effect of y on 
inflation is discussed in section VI. The time at which dark energy fails to get renewed is 
calculated in section IV. B. 



III. CREATION EQUATION 



In this section we describe the first new element of our model: Dark energy implies 
a flow of energy from a background, which in our world may be considered as stim- 
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ulated emission. Stimulated emission scales as the square of the relevant phase space 
number density, and in our model we identify this phase space density in analogy to 
the holographic principle. Figure [T] shows the process of coherent emission of gravitons 
from the background. We picture the background as made up of particles at the high- 
est energy conceivable, energies at which all natural forces combine; we will call these 
all-encompassing particles gravitons to emphasize the dominant nature of gravity. Since 
gravitons are bosons, let us consider stimulated emission and determine the the expres- 
sion for the number density for these bosons. For photons, which are also bosons, the 
basic equation for stimulated emission is the Boltzmann equation describing the scatter- 
ing of photons on electrons, in this application commonly referred to as the Kompaneets 
equation, e.g. [9]. For a discussion of the Kompaneets equation see Appendix B. 




FIG. 1: The stimulated emission of photons 

Whenever a black hole is formed, a large number of gravitons is generated, and we will 
argue below that it is plausible that most of these are produced in the very last collapse 
phase, just prior to making the black hole proper. So a shell of gravitons, which we call 
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the "shell front" , is produced and propagates outwards as shown in Fig|2j 




FIG. 2: The stimulated emission of gravitons in a shell 

In the following sections the Lorentz factors for the shell fronts and for black holes are 
derived, and an expression for the critical density is obtained. These expressions are then 
used to determine the rate at which gravitational energy is created in our universe. 

A. Stimulated emission from the background gravity brane 

We propose a model in which energy is transferred from the strong-gravity (strong- 
gravity) brane (see Eq.(|3])) to our brane via resonance in gravitons, that is by stimulated 
emission. Here we first give the microscopic derivation for all three cases, dark energy, 
black holes, and inflation, before moving on to integrate the microscopic derivation into 
a cosmological context. 

We are using the concept that in a 5D-space the phase space distribution is Planckian, 
with all the momenta in the fifth dimension being maximal, at Planck level. This modifies 
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the 4D phase space distribution such that the momentum in the fifth dimension acts like 
a chemical potential making the distribution essentially flat for all energies up to Planck 
level. Interestingly, this implies that in such a case the Fermi and Bose statistics are very 
nearly the same in 4D-space. 

The Planckian distribution in five dimensions is 



e {hu>/(k B T Pl )} _ I ( 17 ) 

with 

(M 2 = (cp x ) 2 + {cpy) 2 + {c Pz ) 2 + c 4 m 2 Pl , (18) 

where the gravitons have a momentum with components p x , p y , and p z . So under all 
reasonable circumstances 

hU ^1. (19) 



ksTpi 

The energy density of such a distribution is naturally of order m P ic 2 /l Pl . We posit that 
this phase space distribution can never be exceeded. This distribution is modified at 
low energies. On the strong-gravity brane from the Higgs mass on down the distribution 
function is modified due to the rest mass of the Higgs. This lowers the value of the 
distribution function. 

We orient the coordinate system such that any motion or gradient of a potential is 
aligned with z. We associate wavelengths via = h/p x , \ y = h/p y , and \ z = h/p z . We 
define the Lorentz factor T P i with the motion which is required to bring the energy cp z 
up to Planck energy level. 

We then assume that in the proper frame, to be defined below for each maximal 
rate of energy transfer from the gravity brane to our brane is given by the rate gravitons 
in resonance move through the surface of the graviton, multiplied with the fraction of 
maximal momentum phase space. 
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This maximal rate of energy transfer is modulated by the square of the coordinate 
space density on a surface in Planckian units. The transfer is derived from a rate at 
which a graviton on our brane becomes resonant with gravitons on the gravity brane, also 
on a surface. The square of phase space density is obtained when producing a graviton 
each time, since we have an equal number of gravitons on the gravity brane and on our 
brane, all of which are in resonance. 

At the end we reformulate this into some simple rules. 



In this section we describe the shell of gravitons emitted when a black hole forms and 
travels through space. Further below we argue that the uncertainty in space and time 
implies a very large Lorentz factor for this shell front; the speed of the shell front is just 
below the speed of light. This Lorentz factor compares the scale of the universe with the 
Planck scale and is given by 



where H(z) is the Hubble parameter, and r(z,z*) is the spatial integral from redshift z 
to 



In concept this Lorentz factor is similar to the Gibbons-Hawking temperature. 

This means that in the proper frame the ^-component of the momentum is increased 
by a factor of T D , and its length scale is reduced by the same factor. 

We now consider the area of a graviton in the proper frame. In the comoving frame of 
the shell front the rate of resonant gravitons going through can be written as a multipli- 
cation of several factors, the first of which is an emission transfer rate 



B. The shell front 




Pl 



(20) 




(21) 



TT\ X XyCT D 



(22) 
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The speed of light appearing here is similar to a corresponding expression in radiative 
transfer, where emission can be written as energy density times c/(4ir). We will multiply 
with the energy density below. The time step here is, however, the Planck time, and we 
can rewrite the speed of light as Ipi/rpi to obtain 

nX x X y lpiT D — , (23) 
Tpi 

which appears as a volume transfer rate due to the finite time step. The second is the 
fraction of momentum phase space in resonance 



ipi r D 



nXxXy Tpi 

The third is the energy density of the gravity brane 



(24) 



(25) 

L Pl 

Transforming everything to observer frame takes out Y 2 Dl one for energy and one for rate. 
The total transfer rate, which is obtained by multiplying all fo the terms above, gives an 
energy rate per patch of 

mptc 2 1 ,^„ s 
~~ pr — • (26) 

1 PI Tpi 

Multiplying this with the number of patches and the Hubble time yields an estimate of 
the total energy transferred from the gravity brane to our brane per black hole. This 
gives the maximum rate of energy transfer. 

For this to reproduce dark energy as an energy transfer from the gravity brane this 
calculation shows that the extreme energy density on the gravity brane is required. We 
will go through the numbers further below. 

However, this effect is modulated by the proper phase space density squared to obtain 
a graviton, which we derive next. 
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C. Critical phase space density 

The surface area of the front is A, and the number of gravitons is N. In terms of these 
quantities the surface density per Planck area, which is the filling fraction of gravitons, is 

N T D . (27) 

This can be rewritten as, following the line of reasoning above in analogy to radiative 
transfer, 

iV^P — • (28) 
A r Pl 

We recognize that this measures the rate at which resonant gravitons go through the 
surface in the comoving frame per Planck volume. This becomes the critical case when 
this rate reaches unity. For stimulated emission we postulate that this is the critical 
phase space density and for the creation of a pair of gravitons this normalized phase 
space density must be squared. 

This definition of criticality is closely related to the concept of the entropy of black 
holes, where also the surface density is of order unity [T9H28] . 

We obtain a square here due to the fact, that we are using the same number of gravitons 
on our brane as on the gravity brane, which is analogous in the language of Feynman (vol. 
Ill, eq 4.28 on p. 4-8; [29]) of having n photons together with n atoms. In the present case 
Feynman's atoms are analogous to the graviton distribution at momenta at or above the 
resonant graviton under consideration on the gravity brane, and the photons are analogous 
to the gravitons on our brane. In our picture there is full resonance, so all gravitons on the 
shell front involved in the energy transfer are in resonance with gravitons on the gravity 
brane. But much of momentum phase space on the gravity brane extends beyond the 
resonant gravitons, all at the limit of momentum phase space can bear. The extremely 
small momentum phase space factor controls this (of order 10~ 105 ). 
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The number of patches can be written similarly as 

7T \ X Ay 

This expression can be reformulated as the local rule that at most one resonant pair of 
gravitons is created per patch (n\ x \ y ) per Planck time in the proper frame. Transforming 
to the proper frame implies ((n\ x \ y )/T D ) is the size of the patch. As noted above the 
single factor of T D arises from the condition that the argument derives from a rate. 

By similar reasoning the surface density of gravitons is given by a surface density in 
the observer frame multiplied by the Lorentz factor to put it in the proper frame, 
since this similarly derives from a rate. 

When we generalize specifically we will add the required redshift factors. 



D. Black holes 



Black holes have a deep potential well, with a Lorentz factor of 

T BH = ¥™ . (30) 
mpi 

The momentum p z goes up to Planck scale, but the lateral momenta remain unchanged. 

The area for a single graviton is given by n\ x \ y , which in this case corresponds to 
the horizon of the black hole. The final result is that we obtain for the three-momentum 
space fraction, after combining the three coordinates, 

I pi I pi , „ 

-1, (31) 

A X Ay 

which gives an emission rate in the observer frame whose order of magnitude is 

rn^Vmpil 

M BH J T Pl 

One power of the black hole mass comes from the energy of the graviton in the observer 
frame and the other from converting the rate to the observer frame. 
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The surface density of gravitons is quasi-automatically critical for a black hole. The 
entire surface is only one single patch. Thus the same simple rule described above can 
be applied: at maximum one resonant pair of gravitons is created per patch (irX x Xy) per 
Planck time in the proper frame. 



E. Inflation horizon shell 



For the inflation horizon shell the derivation of the Lorentz factor for the shell front 
differs from the derivation of the shell front Lorentz factor in Section 2 in that there is 
no time uncertainty only a spatial uncertainty, because we can not say to better than a 
Planck length, where the horizon is. The details are given further below. In such a case, 
the Lorentz factor of the horizon is given by 

r r infl / qq n 
i- infl — —j , [OO) 

which is a linear dependence. At the end of inflation the value of Ti n fi is about 10 17 ' 7 . 
This implies that in the horizon shell we pull up gravitons from a lower energy in the 
background by this factor. At the end of inflation this is given by 

r ™f l = i n jT~7 — W2 ' ( 34 ) 

fpd 1 + z infl) ' 

On the left hand side we have 



mpic 2 

Fi n fi 



(35) 



and on the right hand side we have 

ch 



l + ZinflF- (36) 
T H C ■ 

These two energies are identical, showing, that on the gravity-brane we reach down to 
the Higgs mass at exactly the same time that we reach the Higgs mass [3UJ [3T] at the 
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observer location in the inflationary bubble. The Higgs mass defines the end of inflation, 
since from the Higgs mass down on the gravity brane the distribution function is modified 
due to the rest mass of the Higgs. This lowers the value of the distribution function, and 
so pulling up gravitons through the Lorentz factor of the horizon shell ceases. The role 
of the Higgs particle in the inflation of the universe is discussed in more detail in section 
VIII.E. 

Here again the surface density on the horizon is maximal anyway, and the patch scale 
is Planckian already. 

Therefore the simple rule formulated above that at most one resonant pair of gravitons 
is created per patch (n\ x \ y ) per Planck time in the proper frame holds here as well, so for 
all three cases. 

F. Rate of Energy Creation 

In the construction of the creation equation we assume that a local comoving frame 
exists in which the energy is transferred from the background. We will show in the 
following section that a Lorentz factor of the comoving frame can be derived from the 
uncertainties in space and time. Now write these rules as a function of redshift, so 
generalizing to the cosmological evolution. 

The creation equation can then be written in the form of energy per unit time created 
in the form of resonant gravitons 

f = ( NGW £k rD ) {^ Td ) EgwTeT ~p' W 
in the comoving frame with proper time r on a spherical surface of radius R s with graviton 
number New (An expression for Nqw as a function of the mass, Mbh, of a black hole 
is given in section IV.B). Here a P i is the Planck area, and r P i is the Planck time (see 
section IX.F.l, List of definitions for the exact expression and the numerical value). T D is 
a Lorentz factor applied to the surface density relative to the observer, and to the size of a 
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patch in the comoving frame, and Te the Lorentz factor as applied to the transformation 
of energy. In the next section we show that for gravitons Td = Te- Eqw is the energy of 
the gravitons on the surface in the observer frame, and A ~ ch/Ecw is the corresponding 
length scale, with the numerical constant needed to make this relation an equality given 
by appropriate matching conditions. 

This macroscopic expression is constructed to be similar to the nonlinear Boltzmann 
equation for stimulated emission in the relativistic case discussed in Appendix A, and also 
the Kompaneets equation describing stimulated emission in the non-relativistic case (see 
[8-10J), but is really based on the microscopic derivation presented above. The first term 
is in the form of an occupation number, and enters the relationship squared. It is the 
density of gravitons per Planck area in the frame of the surface. The second term is the 
number of correlation patches in the surface in the observer frame, and Ts is the Lorentz 
factor which transforms this number to the co- moving frame. This Lorentz factor is 
necessary, since in the shell frame the correlation patch is smaller, and so a Lorentz factor 
>> 1 applies. The third term is the energy in the observer frame, again transformed to 
the co-moving frame a Lorentz factor, and the last term is the rate in the surface frame. 
To transform to the observer frame we take out the Lorentz boost of the graviton energy, 
and decrease the rate by the same factor. 



The expression in Eq|37]will be utilized in the following sections to describe the phe- 
nomena of dark energy, black holes, and inflation. 



IV. DARK ENERGY: SHELL FRONT 



We have argued above and will discuss further in section V that in the last moments 
of black hole formation a shell front, consisting of many gravitons, is produced. To 
describe its motion, we derive a Lorentz factor for the motion of the shell front given 
by a combination of the spatial uncertainty lpi/r(z,z*) with the Planck time uncertainty 
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TpiH(z) relative to the expansion rate of the universe; this is in a direct analogy to (2j |3] 
in which a combination of the scale of the universe with a Planck scale is used. r(z, z*) is 
the distance scale of the universe, z is the current redshift, and z± is the redshift when the 
soliton front was born in the formation of a black hole, or analogously, in the merger of 
two black holes as in Fig|3j We emphasize that we are assuming that most of this energy 



At 



/ BlackHole Merger 

FIG. 3: Propagation of a shell front created by the formation of a black hole or the merger of 
two black holes. 

emerges in the last few characteristic collapse times, so that the shell front has a thickness 
of just a few wavelengths in the observer frame, a characteristic which we will use. There 
is a spatial uncertainty giving one Lorentz factor, and a temporal uncertainty giving 
another Lorentz factor. The combination of the two uncertainties decreases the effective 
Lorentz factor, since for an increase the front would be outside of resonance with the 
background. These gravitons emitted into our world also have a quasi-Planck spectrum, 
with a characteristic temperature - given by the original black hole which creates the front 
- which decreases with the expansion of the universe. 
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A. The Lorentz factor of the shell front 

In the following we use the analogy with the emission of a photon simultaneously with 
the emission of a particle of ultra high energy. After traversing a distance of the scale of 
the universe the time delay of the particle with respect to the photon is 



^ttrav o 2 ^trav (^8) 



1 

2 ifjHECR 

where Jurecr is the Lorentz factor of the particle, t tra v the travel time, and At tra v the 
time delay. We will use the same expression for the delay below. 

First we have to ascertain how precisely we can localize the shell front in the observer 
frame. The shell front has a thickness of several gravitational wavelengths, say 3 £^0.5 

\ q z M bh 1 + 2* 

*gw — As sh Q^Lpi — — . l^yj 

mpi 1 + z 

However, since we have a very large number of gravitons, the precision to which we can 
measure length is in turn given by 

^gw M B h 1 + z* Slpi 1 + z 

V N GW V N GW m Pl 1 + z V4:1T 1 + 

which means that the precision is of order lp[. It can obviously not be less than Ipi, 
the limit is in fact exactly lp\. For our adopted parameters New — 10 97 . In fact, this 
suggests that this equation misses a factor of 3/a/{4~t} ~ 1, which we will use below. 
This factor is similar to an equivalent factor of order unity in optical diffraction theory. 

Any modification of the velocity of the shell front implies a speed below the speed of 
light. The difference between this velocity and the speed of light can be derived using a 
random-walk analysis. There are two uncertainties, one spatial and one temporal. If the 
velocity differs from the speed of light, it can only be lower. For any motion slower than 
the exact speed of light c we have a speed (3c, and the distance traversed towards the light 
front can be written as some distance d. The spatial uncertainty of the distance could be 
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positive or negative, but the distance is just slightly smaller than the distance to the light 
front. 1 — (3 can be written also as 1 — l/(2r 2 ) for T » 1, where T is related to (3 as 
T~ 2 = 1 — (3 2 . Since the velocity is (3c, the advance is (3ct, where r corresponds to the 
age of the shell front, which is some time of order the age of the universe. The relative 
delay in the distance travelled is 

L cr. (41) 



2r 2 

The uncertainty in each step is the Planck length Ipi. This spatial step-size can be either 
added or subtracted, and thus going from one time step to the next is equivalent to a 
one-step random walk. Therefore the length squared d 2 changes on average as d 2 + l 2 Pl , 
where d is the total distance travelled. Due to the one-step random walk, the uncertainty 
in the advance is 



This implies by comparison 



1 I 2 

2d 2 1 ] 



1 l 2 

— = (43) 
2T r , 2 2d 2 ■ [ 6) 



d 

So we can identify the corresponding Lorentz factor as Yd = d/lpi. We can identify d 
with the distance scale, r(z, z*), in our universe since we are considering the propagation 
of a light front (see section IX.F.2, List of definitions for a definition of r(z,z*)). We are 
making a local comparison, so there is no extra factor of (1 + z±). The Lorentz factor in 
this case is 

r„ = . (44) 

In addition there is a corresponding uncertainty in time itself, which can also be written 
as a propagation slower than the exact speed of light. However, time can go in only one 
direction, so we do not have to add squares. The temporal contribution to the Lorentz 
factor is determined from the probability that each step occurs. This contribution is thus 

1 H(z)t pi 



2T 



2 " o ( 45 ) 
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Therefore we obtain 



r t = {H{z)r Pl }- 1 ' 2 . 



(46) 



The final expression is obtained by combining these two Lorentz factors to obtain a 
common Lorentz factor for the shell front. Adding the velocities would boost the Lorentz 
factor to the point where a normal graviton is outside resonance even with a Planck energy 
particle in the background Planck sea. Therefore the velocities should be subtracted, 
which by the usual rules of relativistic velocity transformation is equivalent to dividing 
the larger Lorentz factor by twice the smaller one. The final Lorentz factor is 



This is the Lorentz factor of the shell front in the frame of the freshly born black hole. 
Typical numbers are V ~ 10 31 for the Lorentz factor, which implies that the velocity of 
the shell front is 1 — 1CT 62 - 3 times the speed of light. 

Extrapolating the behavior of the Lorentz factor into the future yields an asymptotic 
constant, as H(z) approaches a constant for (1 + z) — > 0, and so r(z,z+) (see section 
IX. F. 2, List of definitions) also approaches a constant. 

This derivation should be applicable to any signal, that is, to any change of a wave or 
massless particle travelling with the speed of light. A change in the wave or an increase 
or decrease in the flux of particles constitutes a signal and thus corresponds to new 
information. The derivation shows that any such signal always travels at a velocity just 
below the speed of light < c, just like a particle with mass. However, then we have to 
ascertain again, whether we can actually reach the Planck limit precision, and the answer 
will generally be no. Therefore the limit comes down to the wavelength of the energy 
packet divided by the square-root of the number of coherent such packets N; this clearly 
says that this argument works only for Bosons, not for Fermions. This gives for the 
uncertainty, using the matching condition above which requires the extra factor of y/An 



r 



1 r(z,z+) 

2 I pi 



{H(z)r Pl } 



1/2 



(47) 



mm 




(48) 
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where N is the number of Bosons of individual energy E\ which act coherently. The total 
energy is then E% = NEi, and we obtain for the Lorentz factor 

1/2 

r, = I £r ( — ) i • (49) 




This should also hold for sub-atomic scales, but it requires a very large number of 
coherent Bosons to be of any interest. We will explore the consequences which may 
arise from this result elsewhere, but will point out one consequence of interest in the last 
section of the paper, section VIII. For a single wave, this energy corresponds to a mass 
of ~ 10 -84 ' 3 g, or, for all of the waves in one soliton, to a mass ~ 10 97 times higher, 
~ 10 12 ' 6 g, numbers extremely small compared to the mass of the black hole just born, 
which is 10 6 5 M Q = 10 39,8 g using the numbers adopted. 

B. The Shell front 

Using the Lorentz factor derived for the shell front, we can work out the rate of dark 



energy production. We have verified that the occupation number in Eq.(37) (for an exact 
definition of this number see Eq.(51)) is close to unity for our adopted parameters. For 
more massive black holes it would be even larger; for significantly lower mass black holes it 
would be much less than unity, and so these black holes would not contribute significantly 
to dark energy production. 

Inserting the number of gravitons Nqw given by the creation of a black hole or the 



merger of two black holes, into Eq,40 above we obtain (spin J = approximation) 



M BH \ 2 fl + z^ ] 



Now = 4vr — — - . (50) 
V m P i J \1 + z J 

Ideally we would integrate the creation equation, and self-consistently derive the time 
dependent dark energy contribution. Since this is beyond the scope of the present inves- 
tigation, we simply determine the power-law dependencies. The exponent of the redshift 
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term follows from the requirement that we obtain self-consistency at the end. No other 
simple power-law dependence would allow this. This approach should be considered as a 
consistency check, since we do not propose a more general equation, and then derive the 
specific expression. However, we show below that dark energy does end, far in the future 
of the universe, so postulating this expression goes beyond reproducing known observa- 



tions. The first term in the creation equation (Eq. 37) is the "occupation number" N occ 
squared, 

"~=C^r D y (5i) 



V ^d 2 L / 

This number is of order unity (using our adopted parameters it is N occ ~ 10 a5 ) and for 
Td = T it is "critical" as required by our approach. For the definition of the luminosity 
distance di see section IX. F. 2, List of definitions. 

For the Schwarzschild metric the scale of the resonant correlation patch is given by 

1 /7T\i/a G N M BH 1 + z, 

2 V 3 J c 2 1 + z 



* = ;m (52) 



where we will justify the numerical factor (7r/3) 1//2 from a matching condition below. The 
factor of 1/2 is appropriate in the case that we identify the diameter of the Planck shell 
with a half- wavelength, that is using adjacent nodes of a wave. A numerical value is 
A ~ 10 13 cm, again using our adopted parameters. The number of correlation patches on 
the shell front in the comoving frame is about 10 66 4 , using our adopted parameter values. 
The resonant energy of a graviton is 



2 mpi 1 + z 



E G w = m Pl c —— — . (53) 

8tt M B h 1 + z, 

In the comoving frame (i.e. after multiplying with V) this is about 10 +1 erg. Finally, 
using the expression for the Lorentz factor from above and taking out the factor V 2 to 



move to the observer frame, Eq.37 becomes, 

cIEt 3 , , 2 ( 1 + z+ 



dt 2 \l + z 



3 



M BH c 2 — — H(z). (54) 
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Inserting all the expressions then gives the redshift dependence of dark energy, once folded 
with the density of black holes, and their redshift dependence. 

We can work out directly the required rate of energy production, and then match the 



result to Eq.54 to obtain the values of the numerical factors in our approximations. What 
we need per black hole is given by 

j t {N GW E GW ), (55) 

which is 

3 , r 2 ( 1 + z,^ " 



2 M BH c 2 ( ' H ^ z) - (56) 



Matching the expressions in Eqs. 55 and 56, and using the expression for \gw m Eq.52 
the factor is found to be ~ 1.02. Setting NqwEgw ~ Pde, we see that the 

energy source in our universe scales with PdeH(z). We will use this below when we will 
employ the energy-momentum tensor and the Einstein equations. 

It is straight-forward to verify that a shell front cannot just collapse into mini-black 
holes. Such a front does not have the critical density to do this during its formation and 
never reaches the critical density. As it expands its density in its own frame scales with 
r _1 , since the scale of the radius enters the Lorentz factor. At the same time its thickness 
scales also with So the two time scales never cross. 

Combining this with the observed space density of supermassive black holes today 
cancels all redshift dependence and gives a numerical value for dark energy, which is 
consistent with observation to within the rather large error bars, as shown in the next 
subsection. 

For a more detailed discussion of the conditions at the end of inflation see Appendix 
C and D. 
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C. Observations 

The relations derived in the previous section for the rate of dark energy production 
can now be used to obtain a value for the dark energy density, which can be prepared to 
the observed result. For simplicity we use the same values of the parameters as in Section 
IV. A. However, we initially keep these parameters arbitrary in the expressions in order to 
show the parameter dependencies. 

When a black hole is created or two black holes merge, the maximum efficiency of 
converting rest-mass into energy is 50 percent. This efficiency is reached when either two 
black holes of opposite, maximal spin and of equal mass merge [32], or when a black 
hole forms out of a collapse into a non-rotating black hole [S3]- We use this limit here 
and address the cases of rotating black holes, orbital spin and mergers of black holes of 
unequal mass in a future investigation. 

The total energy is the number of gravitons per black hole times the energy of each 
graviton times the density of black holes 

„ fM BH \ 2 fl + zA 4 1 fl + z m 2 Pl c 2 \ M , 3 

V m Pl J V 1 + Z J 87r V 1 + Z* M BH J 

To obtain an estimate of the dark energy density we set Mbh = 3 ■ 10 6 M & and z* = 50; 
The density of black holes at creation, using today's value, is tibh,o = 10~ 2 2 Mpc~ 3 [34] . 
The original density of super-massive black holes could be substantially higher, e.g. by 
a factor of order 8, if black holes grow mostly by merging (e.g. [351 [36]). In addition, 
the observational uncertainty in the number itself is known to be a factor of 2.5 (one 
sigma), see |34J. The black hole mass could perhaps be 1O 7 M , giving another factor 
of 3. Furthermore, the creation redshift z± could be larger or lower than our adopted 
example of 50. Using the estimate in [37] of the first stars forming at redshift 80 at 
most, super- massive stars seem unlikely to form earlier than at a redshift of 70, since the 
difference corresponds to a few million years, which is the life-time of very massive stars. 
This would modify the number by a factor of 4. In view of all these uncertainties the 
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parameters allow the match with dark energy of 10~ 8 ' erg /cm 3 ; using today's parameters 
naively without any correction yields a value for the density two orders of magnitude 
lower; this also matches what is usually called the entropy of the universe. Since the 
largest correction here derives from assuming that super-massive black holes grow by 
merging with each other, this seems to be an implied consequence. We do not use the 
black hole density in the creation equation; we follow the evolution of the shell front of 
a single black hole. In order to obtain the observed vacuum energy density today the 
original number of black holes would have to be ubho = 10-°- 1 M^ 6 . 5 Mpc- 3 , scaled to 
today. This value assumes that all other parameters are unchanged. In the total energy 
density the dark energy contribution scales with the total mass in super-massive black 
holes. Dark energy has grown by a factor of about (1 + z+) 3 ~ 10 5 since it was first 
created. 

One check here is that the efficiency of turning normal matter into the first 
super-massive black holes should be small [38J. We use the numbers of ubh,o — 
1CT ' 1 M^hg 5 Mpc~ 3 as the original black hole density scaled to today, and consider a 
black hole of 3 • 10 6 M Q made out of baryonic matter. This estimate uses for comparison 
the entire baryonic mass out to the next sphere of influence of another similar black hole. 
At the adopted formation redshift of = 50 this sphere has a radial scale of about 
lO^l + z^Mpc = 21kpc. Since the specific volume scales linearly with the mass of 
the black hole, dividing the black hole mass by the specific volume renders the efficiency 
independent of the mass. Furthermore, the efficiency is also independent of the formation 
redshift adopted, since we use numbers referring to today. This gives an efficiency of 
3 lCr 5±0 ' 4 , which, considering the original observational uncertainties [31], matches the 
extremely low efficiency found elsewhere. 
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V. BLACK HOLES: PLANCK SHELL 



When a black hole is created (e.g. [33 } l39H4"T] ) . or when two black holes merge [32], we 
argue that all of the rest mass energy can go into the black hole and partially emerge as 
gravitational waves, neutrinos or electromagnetic emission. The limiting case is that half 
of the energy goes into gravitons. This limit is valid for a Schwarzschild metric, and also for 
two equal mass black holes of maximal and opposite spin in the plane of rotation merging 
(see also |42j). Although our metric is not Schwarzschild, we nevertheless assume the same 
division of energy into gravitational waves and neutrinos or electromagnetic emission. 

This can be derived from our creation equation Eq|37} Consider the collapse of a 
massive configuration towards a black hole; this implies that we need to include the 
Lorentz factor corresponding to the potential depth, and further below we will argue that 
this is Mbh /rripi- This necessarily corresponds to a non- vacuum solution. We also want 
to get the total amount of energy coming out in the burst, and so we integrate in the 
observer frame by multiplying with Tp\. We note that at the moment when the black hole 
is born, the first term in the creation equation is unity times T 2 D , the second term is unity 
times 4 To, the the third term is the energy of a gravitational wave in the observer frame, 
Egw ~ times Ye- All this is in the moving frame of the collapsing black holes. 

Here as seen from an observer =1, but Te = Mbh/t^pi- The total energy coming 
out (over Tpi x Te) is 



however, for convenience we assume in the following that the numerical factor in our 
simple approximation is just 1/2, and not unity. This supports the notion, that at the 
formation of a black hole a large number of gravitons are emitted, matching in total of 
order half the rest mass energy of the black hole mass, seen from a distant observer. 

However, this entire argument is more of a consistency check again, since the actual 
value of the emitted energy is obtained by integrating over the progression of the Planck 



Ebh 



^ M BH c 2 , 
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shell, as it first forms, then grows in the collapse of a massive object, see [13]. At each step 
when the Planck shell has enlarged by a Planck area in its own frame, it can emit a graviton 
per Planck time. This is rather similar to the initial formation and subsequent growth of 
a shock wave in non-stationary gas dynamics, which also separates causal regimes, pR) . 

We will deal with black hole spin elsewhere. As a guiding principle in the following we 
adopt this limiting efficiency for the formation of the first black holes, and assume all of 
the radiated energy to be available for gravitons. At very high particle energies we take the 
graviton to incorporate all natural forces. The implication is that at the thermodynamic 
limit gravitons are generated. We identify what is normally called "entropy" of black 
holes with gravitons, both far from the black hole, and also just outside what is normally 
referred to as the horizon. It is this "shell front" of gravitons propagating through the 
universe from each black hole formed (or from each black hole merger) which provides the 
connection of black holes to the cosmos. However, the energy of the gravitons depends 
on the black hole mass, so the energy of these gravitons will generally be different for 
different black holes. This is also different from current approaches (e.g. [7]), where these 
two populations of gravitons do not exist. 

A. Planck shell 

We propose that a black hole can be pictured as a "Planck shell" full of gravitons 
at potential depth MBul m Pi i n terms of a Lorentz factor. This is rather similar to the 
concept of a "stretched horizon" [271 EH] • This shell is an impenetrable barrier, and the 
ensemble of these gravitons constitutes the black hole as shown in FigJIJ Space-time ends 
at the inner surface of the shell; there is no physical meaning to space-time coordinates 
in the region enclosed by the shell. This concept was also later suggested in terms of a 
"firewall" by Almheiri et al. |46j . 

These gravitons have a quasi-Planck spectrum. The energy density on the Planck shell, 
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Black Hole Shell 



FIG. 4: The Planck shell for a black hole 

Ebh, is given by the Planck limit multiplied by the ratio of the Planck mass to the black 
hole mass 

mpi mpic 2 

£bh = — —^3- (59) 
M BH Anl Pl 

The creation equation above requires that energy be pulled from the background, so the 
black hole does not evaporate but very slowly increases in mass. 

The derivation of the Lorentz factor in terms of the black hole mass begins with the 
gravitational potential. The invariant line element for a curved space-time background 
can be expressed in terms of the gravitational potential, $(f), as 

rfs 2 = _ e 2*(f)/c« c 2 dt 2 + g ,, dx i dx j 



The ratio of the frequency of radiation emitted at a point, r~s, at the Planck shell surface 
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to that at a distant observation point, r~o, is given by 

vo y/(9oo(rs)) 



(60) 



The Lorentz factor V BH is related to the frequency ratio in eq. [60] by the expression 

— = v BH (6i) 

and Ybh is related to the goo( r ) component of the metric tensor by the relation given 
above. The form of T B h which gives consistency can be written as 



BH 



f(r) (r - r s ) (62) 



where f(r) is defined by the requirement that T BB has the proper values at infinity and 
at the horizon. We use the ansatz 

f(r) = ^ . (63) 
r 

This results in 

, r s + lpi , . 

fo = —j (64) 

The ^oo element of our metric tensor element then has the form 

rs ^ 2 



and the invariant line element is 



9oo = (I- ^ ) , (65) 



ds 2 = -(1 - ^) 2 dt 2 + - l -— dr 2 + r 2 dO 2 + r 2 sin 2 (^) d<\> 2 . (66) 
2r (!-27) 2 

This results in a Lorentz factor of the required form 

r s M B h 

L BH — 7zy~ — • 

l ipi rripi 
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This metric is of the extreme Reissner-Nordstrom type with the charge Q obeying 

G N Q 2 = (G N M/c 2 ) 2 . (67) 



We are not proposing the metric in Eq. 66 as the correct form for the metric at a 



Planck surface. We want to show only that a Lorentz factor of the required form can 
be obtained for non-vacuum solutions of the Einstein equations. The actual form of the 
metric near the Planck surface is of course determined by the interaction of the surface 
with the background energy sea and with any other fields which may be present. The 
determination of the true metric is beyond the scope of the present investigation but will 
be the subject of a future investigation. 

The filling condition of the occupation number for the Planck shell is obviously satisfied 
(Td = 1 and Te = Mbh/^pi)- In resonance the graviton on the Planck shell is 
exactly the Planck energy in the Planck sea, the background. We also note, that laterally 
the graviton has the size of the circumference of the black hole horizon, and radially 
the thickness is the Planck length. A similar model with an extreme anisotropy was 
investigated by t'Hooft, [17] . If all of the gravitons on the Planck shell are in a coherent 
state, i.e. all of the gravitons are aligned, our model reproduces the spin of a Kerr black 
hole. 

Using the creation equation, the occupation number can be shown to be of order unity 
independent of our standard parameter choices. Creating one Planck graviton per Planck 
time on the Planck shell implies that in the observer frame the luminosity is locally 

1 mpic 2 ( rripi \ 2 , . 
L = — — — . 68 

2 t px \M BH ) 1 ; 

The factor of 1/2 is an approximation as explained in section IV. D. Half the energy 
obtained from the background field is radiated away and half is absorbed, increasing the 
mass of the black hole. There is no net evaporation of black holes in this approach, 
although the Hawking process still happens [HH5T]. However, the Hawking process does 
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not drive the evolution of the black hole, it is actually much weaker than the process 
described here, the interaction with the background. This implies that primordial black 
holes all grow to the point where their growth time scale becomes the age of the universe. 
This allows rather strong limits on the initial number of primordial black holes [52H51] . 
as we will discuss further below. 

The typical energy of the graviton coming out in the observer frame is given by 



To a distant observer this is measured as a redshifted Planck spectrum. 

A black hole forms, when a Planck shell is created at some finite radius in a collapsing 
configuration. Therefore we have to question whether at that point we introduce a new 
variable in terms of how many gravitons get generated as a function of the enclosed 
mass at the moment at which the Planck shell is created. However, since in the end the 
black hole cannot have any extra parameters, such as this specific radius, the graviton 
generation must be independent of any extra parameters. A fortiori graviton generation 
must be the same as if the black hole were to start growing from a Planck mass. The first 
formation of a Planck shell, and its subsequent growth are just a function of its mass. 



The rate at which gravitons are emitted from a black hole in an expanding universe 
without dark energy can be calculated in a semiclassical approximation. Following Parikh 
and Wilczek [55] we calculate the emission rate for a graviton of frequency u by calculating 
the action in the WKB approximation for a black hole in an expanding universe. The 
graviton emission probability for a black hole in an expanding universe can be calculated 
from the McVittie metric [SS] 




(69) 



B. Graviton emission rate in an expanding universe 
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where 



m . (7i) 



fi(t) 2' 

The semiclassical emission rate is thus proportional to (for a more detailed derivation of 
this see Appendix E) 

r em ~e" 235 = e - 4 ™WM 5 "W 2 ) (72) 



Neglecting the term quadratic in u, the exponential in Eq.(72) becomes a Boltzmann-like 



factor with a temperature which decreases with increasing time. This is consistent with 



our model of ever-cooling black holes. The rate expression in Eq.(72) does not take into 
account the cooling due to the interaction of the Planck surface with the background. 
This interaction creates further growth, i.e. cooling, of a black hole. 



C. Accreting black holes 

When black holes accrete, they do not pull any additional energy from the background; 
however, the accreting matter produces a mix of electromagnetic emission, neutrino emis- 
sion, and graviton emission in addition to allowing the black hole to grow. The graviton 
component of the emission is usually dominant, and produces a general gravitational wave 
background, which may be just below the current upper limits [57]. This additional grav- 
itational wave emission, which is some factor of order < 10 above the electromagnetic 
emission, provides a universal background of gravitational waves. 

We can vary our creation equation by assuming that a particle, say a proton, accretes 
and gains energy by the potential depth, thus gaining a factor of Mbh/^pi (using zero 
redshift for simplicity). The energy will rapidly exceed the Planck energy, creating addi- 
tional Planck energy gravitons. We reiterate that for a black hole, both the occupation 
number and the number of correlation patches are unity. This implies that the output 
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becomes (e is the fraction of the accreted energy which is re-radiated) 

dE f mpi c 2 mpi \ ( M bh t pi M bh 
dt V t pi m bh J I rn Pl m P/ 

where the first term arises from the normal production of gravitons, and the second term 
arises from the infall of matter, leading to the splitting of the energy into Planck level 
gravitons, and the resulting multiplicity. We choose e to be one-half to describe the 
outgoing graviton emission in our simple approximation. This equation then reduces to 
the trivial relationship 

dE ._ . 

M = £Mb " C (74) 
Baryon number, charge and lepton number are conserved in the creation of new particles, 
but baryons and leptons constitute only a minute fraction of all the emitted and absorbed 
particles. 

In [31] the energy density of such a background was estimated. The mass of black 
holes contributing to the energy density of the universe can be written as 

n BH = 2 • l0- 6±0 - 40 ( 75) 

Assuming that energy is emitted only during accretion, and that growth of super-massive 
black holes is mostly via this accretion [58] 159] . the energy input into intergalactic space 
is given by 6 ■ io~ 15±0 - 40 eo.3 erg cm~ 3 , using an accretion efficiency of 0.3 to produce 
gravitational waves. For a Schwarzschild metric the maximal number is 0.5, with about 
0.06 going into electromagnetic emission. So the maximum fraction of the energy which 
can appear as gravitational waves is 0.44. For a Kerr metric the electromagnetic part 
could reach 0.43, so the fraction of energy appearing as gravitational waves could be very 
much less [60] . 

In terms of critical energy density the density in Eqj75] transforms to 

n GW , S BH = 6 • 10- 7±0 - 40 e .3 f eff (1 + z acc )~ A (76) 
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with the index GW, SBH signifying the gravitational wave emission due to accretion by 
super-massive black holes. / e // is the fraction of the final mass that is due to accretion 
after formation, and z acc is the redshift for the epoch, when most black hole accretion 
occurred. If black hole growth is mostly due to accretion, the f e ff will be very close to 
unity; if black hole growth is mostly by black hole mergers, then obviously f e ff « 1 
is possible. From cosmic evolution studies, the greatest merger activity is known to be 
somewhere between redshift 1.5 and 2 approximately, so the redshift term may contribute 
another order of magnitude. There will also be a corresponding contribution due to 
accretion by stellar mass black holes. 

In the following sections we note again that black holes probably grow by merging 
with each other, and in that case the gravitational wave background estimated here may 
be an order of magnitude lower. This gravitational wave background is just the original 
gravitational wave production; it should not be confused with the dark energy background, 
whose energy density is of order 10 5 higher than that of the background estimated in this 
section. 

Since we are not proposing a specific model for the production of this energy density as 
a function of black hole mass and cosmic epoch, we refrain from proposing a spectrum for 
the associated gravitational waves. However a reasonable assumption is that the energy 
density spectrum reflects the mass function of super-massive black holes. Thus it may 
be quite low at the frequencies associated with 10 9 M & black holes, and peak at the 
frequencies associated with 10 6 5 M black holes, shifted by the redshift dependence of 
the accretion history. This specific number of 10 6 5 M is uncertain by possibly a factor 
of 3, as already noted. 
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VI. INFLATION: HORIZON SHELL 

The model for graviton creation described in the preceding sections may also explain 
the inflationary period [6TH63] . During inflation new gravitons are produced at the horizon 
in a "horizon shell" , provided that the newly created gravitons lose their energy adiabat- 
ically with a time scale proportional to the Planck time, independent of redshift within 
the inflationary period. Exponential behavior is common in any approach to inflation. 
However, in our model the structure of the universe during its inflationary period is the 
one seen by an observer. The inflationary period stops when the new gravitons pulled 
from the background reach the particle energy where the background fails to provide 
enough gravitons. In this concept the universe starts with one Planck-energy graviton in 
one Planck volume. This kind of speculation has been explored for many years, e.g. by 
Lemaitre [12] . 

The concept is that near the horizon we also have a shell, which is at maximum energy 
and particle density. Thus it satisfies the filling condition. 

Within the scenario above the creation of new gravitons always occurs inside our 
universe. Each new graviton which is produced stretches the wavelength of the gravitons 
produced prior to the new graviton. This behavior can be diagrammed as 

If the universe has a size given by 

A< = N n lpi (77) 

where 

N n = N R(t n ). (78) 

where Nq = 1 in our model, and this equation represents a 1-D slice through the early 
universe. 

From the diagram the expression for R(t) is seen to be 
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FIG. 5: The early growth of the universe starting from a single Planck cell. Each cell stimulates 
an additional cell of the same energy. The wavelengths of the older gravitons (cells) are stretched 
by the addition of new gravitons. 

R{t n ) = N 2 n n = ^ (79) 

Tpi 

and 

R(t n ) =N e nla2 = N e {tn/Tpi)ln2 (80) 
Thus the initial wavelength 

Ao = l p i (81) 
will increase as more and more gravitons are produced, so that 
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A n = R(t n )\ (82) 

This of course is an approximation because at some point the wavelength of each gravi- 
ton becomes too large and the frequency becomes to small for the graviton to stimulate 
emission of new gravitons at the peak energy of mpic 2 . The older gravitons generate new 
gravitons at the lower frequency and therefore with lower luminosity. 

This is then directly the normal exponential growth 



r = r e t/T "^, (83) 

where we identify tq with the Planck length Ipi, and Ti n fi with the Planck time Tpi to 
within a factor of order unity. As shown earlier for the shell front, this then gives a 
Lorentz factor of 

IW, = {r/lpi} = e*/^-/.}. (84) 

Therefore the apparent velocity of the horizon increases exponentially, a velocity which 
we translate into a Lorentz factor. A shell front with this Lorentz factor, interacting 
with the background will resonate in the background with particles with lower energy, 
and therefore at lower phase space density. The contributions to the counting of the 
powers of the Lorentz factor include the horizon shell, the surface density squared, the 
number of correlation patches, the particle energy, the rate of particle creation, and the 
transformation to the observer frame. In the end we obtain exactly the three powers 
which we require to push particles to the peak of the Planck spectrum in the horizon shell 
by the Lorentz transformation. For the duration of the inflationary period the number 
density and energy density in the horizon shell are exactly at the background limit, and 
stay there until the refilling begins to fail. This compensates the power-law term of the 
Planck spectrum of (fkj) 3 ^-, after taking into account the factor of hu the denominator. 
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In the exact formulation of the background spectrum the exponent of the exponential 
term in the denominator has to be written as (see section II) 



with y time-dependent, and x the normalized graviton energy, defined above in Eq. [TT 
At the end of inflation the dimensionless parameter x has the value 



x = ~ 10~ 18 (86) 

mpic 2 



The dimensionless parameter y is 

y ~ ^-A\ Tpi . (87) 
mpic 

The parameter 2y reaches the same order of magnitude as x 2 at the end of inflation 
provided that k ~ 10 +50 erg/cm or in terms of Planck units k — empic 2 /Ipi, where 
e ~ 1 now in dimensionless units. When y reaches the same order of magnitude as x, 
the slope of the graviton energy spectrum steepens, which effectively brings inflation to 
an end. Gravitons continue to be created, but at much lower energies and in much lower 
numbers. The value of the parameter e ~ 1 can be identified as the point where a phase 
transition occurs; the Higgs boson begins to appear, and the force of gravity begins to 
diverge from the other fundamental forces and inflation ceases. 

The model which we propose for the inflationary phase of the universe is that on the 
horizon shell the energy and particle density are maximal (Planck level), but particles 
lose energy with the Planck time as they propagate away from the horizon shell. 

Going backward in time from today, we know that the dominant energy density ob- 
served today, which unlike dark energy was not created by drawing from a background 
field, should correspond to the maximal particle energy and particle energy density at the 
value of the redshift at which inflation stopped. There are three possible energy densities 
today, which may relate to this epoch: a) The photon density, at around 10~ 12 erg cm~ 3 
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today, b) the matter density, and c) a possible unknown gravitational wave energy density 
at a level below dark energy as we observe it. The energy density at the end of infla- 
tion near the horizon shell is mpic 2 lp 3 = 10 115 3 erg cm~ 3 in our model. To reach this 
level starting from today's dark energy implies a redshift of 10~ 8 (1 + Zj n //) 4 = 10 115 ' 3 in 
units of ergcm~ 3 , which implies for radiation (1 + Zi n fi) — 10 30 . Any gravitational wave 
background below dark energy is so highly dependent on the model, that it is hard to 
pin down. In our specific model it would be the redshifted Planck spectrum, which is a 
very steep spectrum in energy or frequency up to an energy of mpic 2 (l + Zinfi)" 1 . This is 
a frequency of approximately 10 +10 - 5 Hz, at an energy density of about 10~ 12 ' 7 erg cm~ 3 
(using z in fi ~ 10 32 ), corresponding to a contribution to closure density of about 10~ 5 ' 7 
in dimensionless units. This is consistent with canonical models, but in our approach of 
course the inflationary period is completely dominated by gravitons, at maximal particle 
energy and maximal energy density. 

This then yields a redshift since the end of inflation of about 10 30 , crudely consistent 
with other models, which suggest 10 29 (e.g. [641165] ). This simplification also reduces the 
left-over residual graviton emission spectrum. 

Using 10 29 as the value of the redshift at the end of inflation, we obtain for the horizon 



at that redshift a distance scale of 10 28 ' 2-43 ' 5 cm = lCT 15 ' 3 cm for r in Eq 83 Since we 
start from a Planck length, ~ 10~ 33 cm, this implies an increase of the radius through 
the inflationary period by a factor of 10 17 ' 7 , corresponding to about 41 e-folds, and an 
increase in the number of gravitons by a factor of 10 35 4 . Since this is an estimate, the 
error bars are large. 

This corresponds to a horizon shell Lorentz factor of 10 17 ' 7 (see above) at the end 
of inflation. This in turn implies that particles are created from the background with 
an energy at a factor below the Planck energy, that is an energy of about 10~ L4 erg ~ 
10 10 4 eV. This is an interesting result because this is an energy which is typical of the 
energy of particles observed in high energy cosmic rays (e.g. [66]); we note below that 
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this corresponds crudely to the Higgs mass. We emphasize that this is a particle energy 
in an environment at maximal energy density, thus this result has no obvious relevance 
for our environment of today, since the energy density is lower by a factor of more than 
10 100 . 



This concludes the use of Eq.37 to describe the cosmological phenomena of dark en- 
ergy, black holes, and inflation. In the following sections we analyze the predictions and 
consequences of the model. 



VII. PREDICTIONS OF THE MODEL 



Any model of physical processes is of real value only if makes predictions which are 
observable. Does the model described above have any consequences which are testable by 
observations that can be expected to be performed within the next few years? Fortunately, 
the answer is yes. 

In this section we discuss the various contributions to the gravitational wave back- 
ground and the experimental techniques which can be used to detect those gravitational 
waves which constitute the effect known as 'dark energy'. 

There are three experimental tests which can be performed to determine the validity of 
these ideas: i) pulsar timing arrays (PTAs) to detect the gravitational wave background 
which makes up dark energy, ii) Lunar laser ranging to detect the single solitons coming 
through, and iii) ultra high precision clocks to detect the noise due to the individual 
solitons coming through. 



A. Primary source of dark energy - stellar-mass black holes vs. super-massive 

black holes 

Since black holes of any size can stimulate graviton emission from the Planck sea, a 
comparison of the relative contributions to dark energy from stellar-mass black holes and 
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super-massive black holes is in order. A useful piece of information is, whether or not 
the super-massive black holes have a characteristic low mass, and observations strongly 
suggest that they do [MJ I67H70] . So it is plausible to suggest that super-massive black 
holes may have all started with a mass near this low cut-off. 

One question relevant for the entire approach described in this paper is when the 
earliest super-massive black holes might have formed. In |37j it was shown that for 
one specific model of dark matter star formation could be initiated very early (asee also 
[7T| 172])); obviously there may be other models also allowing rather early formation of 
massive stars. Since massive stars usually form in dense groups, stellar agglomeration is 
possible ([731 HI])- It was shown in reference [75] that the correlations between galaxy 
parameters and super-massive black hole mass apparently merge into a correlation with 
nuclear star clusters, with a different slope [761480] . On the other hand, massive stars 
at near-normal chemical abundances have powerful winds limiting their mass to several 
hundred Solar masses ([HI])- For a heavy element abundance close to zero, such winds do 
not exist, and the first generation of super-massive stars can form and produce the first 
super-massive black holes. At about 10 6 M such stars suffer from an instability due to a 
combination of high radiative pressure, and subtle effects of General Relativity ([82, 83J), 
exploding readily. They could then give rise to the first super-massive black holes, with 
masses between 10 6 and 10 7 M , allowing for some larger mass due to infall. Therefore, 
we adopted above as our reference parameters 3 • 10 6 M Q for the black hole mass and a 
redshift of creation of order 50, the value for the Lorentz factor derived above results in 
an "occupation number" , of order unity, as required by our earlier condition. 

This result implies that stellar mass black holes today violate the occupation number 
condition, since they drop out at a very small luminosity distance, of order 10 -12 of 
the luminosity distance today versus redshift of order 50. Therefore the stellar mass 
black holes do not contribute significantly to the steady creation of gravitons drawing the 
energy from the background; on the other hand, they do keep on creating new gravitons in 
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accordance with the creation equation, just at a very low level. However, the gravitational 
waves generated at their creation should produce a background at the corresponding 
frequency, which may be detectable with normal gravitational wave detectors; this is the 
long expected background. 

The above analysis shows that stellar mass black holes do not make a significant 
contribution to dark energy. However, they do produce a gravitational wave background, 
and they do connect very weakly with the background Planck sea. 

B. Gravitational wave background: primordial background 

There is a primordial gravitational wave background at very high frequency left-over 
from the inflationary period (e.g. |84J). The peak frequency given by the horizon shell 
should be near the Planck energy divided by 10 32 , giving a frequency of 10 10 5 Hz, cor- 
responding to the last Planck energy gravitons produced in the horizon shell at the end 
of inflation. However, the spectrum should extend steeply to lower frequencies, with the 
very lowest frequency being the Planck energy divided by the redshift factor as measured 
from the beginning of the universe of ~ 10 47 7 , giving a frequency of ~ 10~ 4 2 Hz. For 
simplicity we have used a delta-function approximation for the background Planck sea 
of gravitons, but in reality the gravitons created from the background have a Planck 
spectrum. There is an extension to even lower frequencies, again with a steep spectrum. 
The energy density of this background is not far below the microwave background energy 
density. 

C. Gravitational wave background: super-massive black holes 

The gravitational wave background whose existence we propose here is distinct from 
the primordial wave background. The gravitational waves in our model arise from the 
production of black holes in the early universe at a redshift of (order) 50. These black-hole- 
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produced gravitational waves maintain a constant dark energy density by continuously 
interacting with the Planck sea background. In the simplistic framework that all black 
holes were created at the same time, and with the same mass, this implies that the 
gravitational wave background peaks near 

w~ = i \ m ^ 2 w 1 - ttt" * 10 " 4 ' 5 Hz M *™* tt~ > (ss) 

o7T II MbH J- + Z* 1 + z * 

with a steep spectrum at lower frequencies, matching the Planck spectrum background, 
and a sharp cutoff at higher frequencies. In Fig. [6] we show a graph giving current limits 
on the gravitational wave background from pulsar timing ([851 EH]). 

In this graph we note that using lower frequencies or higher initial black hole masses 
would approach the existing limits and thus can be ruled out today. We propose that 
extending the observational sensitivity right at the point of inflection of the observational 
limits, near Q ~ 10~ 6 and / ~ 10~ 6 ' 2 offers the best chance to confirm or refute the 
proposal made here. A further improvement of the observational sensitivity by one or two 
orders of magnitude may be required to carry out this check. 

Considering the many different contributions from merger events giving lower frequen- 
cies and smaller redshifts, the low frequency side of the spectrum might be a bit shallower 
than suggested by just the first formation of super-massive black holes. Even matching 
the known dark energy contribution to the energy density of the universe, this background 
cannot be detected today. However the sensitivity needed to detect this background is 
close to the sensitivity of the planned extension of the pulsar timing arrays (e.g. [87,, ,88]). 
With today's limits from pulsar timing arrays it may already be possible to constrain the 
mass of the first generation of super-massive black holes. We predict that this background 
will be detected. 
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FIG. 6: Gravitational wave background limit from pulsar timing (dashed lines: |85l 86j), and 
our inferred gravitational wave background from stimulated emission of gravitational waves from 
the background Planck sea constituting dark energy, using the frequency scaling given in Eq. [88} 
We use as the ordinate the fraction of closure energy density f2, and as abscissa the frequency 
of the gravitational waves /. 

D. Gravitational wave background: accreting black holes 

Accreting and merging black holes emit a large fraction of the accreted mass/energy 
into gravitational waves, creating a gravitational wave background (e.g. [89]). This 
background folds the temporal evolution of the growth of black holes with their mass 
distribution. The contribution of such gravitational waves to the background can be 
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estimated by the current mass density of black holes which have a typical redshift of 
growth. The current mass/energy density is about 10 -16 ' 9 erg/cm 3 . This value has been 
corrected for the redshift assumed above for the creation of the first generation of super- 
massive black holes, with a large error bar, as noted above. There is an additional 
uncertainty due to the lack of knowledge about which process is dominant in the growth 
of black holes - accretion or mergers. In this paragraph we confine our discussion to 
accretion. If the accretion process is less important for growth than mergers, then this 
background is weaker than that of the other processes, possibly much weaker. In addition, 
since most of the growth of black holes was at redshift of approximately 2 or larger, there is 
a downward correction of order 10 - 20, giving an energy density of 10~ 17 ' 9 erg /cm 3 . This 
translates into a fraction of the critical density of about 10 -9 ' 9 , clearly below current limits. 
The proper frequencies correspond to the mass range of super-massive black holes, which 
range from about 1O 6 M to 3 • 1O 9 M . Improvement in pulsar timing techniques may 
allow detection of this background some time within the next few years [87J. In contrast 
to the contribution to the background from dark energy, this accretion contribution to the 
background is steady. The dark energy contribution is divisible into single soliton waves, 
or shell fronts of gravitational waves, which vary in time. 



E. Gravitational wave background: stellar-mass black holes 

The creation of stellar-mass black holes, which are expected to be more common in 
the early, low heavy-element universe [90] , also provides a gravitational wave background. 
Their typical frequency is 

fGW,stars = ^ \ ^ ^ - 10 +2 '° Hz M^, ^ . (89) 
o7T il NLbh 1 + 

The level of this background at the time of formation may have been an order or two 
orders of magnitude above that generated by the formation of super-massive black holes. 
However, stellar black holes do not give rise to the shell fronts that keep up the filling con- 
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dition which is necessary in order to draw energy from the background. The gravitational 
waves generated during the creation of stellar-mass black holes expand just the same as 
any relativistic gas, that is with a redshift factor of (1 + z) A from creation. Assuming as 
a reference an initial redshift in the range of 50 to 80 [37J, their energy density is diluted 
today by about 10 7 . The existing observational constraints ^77\ give an upper limit of 
about 10~ 5 of closure density, so at this limit the mass density of stellar black holes could 
exceed that of super-massive black holes by a factor of 100, and still be allowed. However, 
this demonstrates that a detection may be possible within the foreseeable future with 
laser interferometers. This specific background of gravitational waves should have a time 
dependence just like the corresponding background from the formation or the merger 
of super-massive black holes. In addition, super-massive black holes should provide a 
contribution to a steady background of gravitational waves due to accretion episodes. 



The increase in dark energy is due to the large number of individual soliton fronts 
emanating from the creation and mergers of black holes in the early universe. Each 
individual soliton front carries a large energy, taking into account the original rest mass 
energy of the black hole, and the large subsequent increase in energy due to the creation 
of resonant gravitons by the solitons interacting with the background Planck sea. This 
energy E s can be estimated by 



For Mbh = 10 6 ' 5 M & and z± = 50 this energy is today (z = 0) 10 65,6 erg. This energy is 
distributed over a volume given by Aird 2 L = 10 61 66 cm 2 ; the distance scale corresponding 
to the age of the universe is 10 28 11 cm, and the corresponding volume is 10 89 3 cm 3 . The 
energy density per soliton of is thus 10 -23,7 erg cm~ 3 . The overall change in dark energy, 
if it were to just relax without any re-creation, would correspond to 10~ 8 ergcm~ 3 4i7(,2) 



F. Lunar laser ranging 




(90) 
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or about 10 _25 ergcm _3 s _1 . With these simplistic numbers we obtain a front coming 
through every 20 seconds, but the transient time is of the order of 1000 seconds. Therefore, 
any given moment there are roughly 50 uncorrelated shells going through. This time 
between shell fronts scales linearly with Mbh and also linearly with (1 + z*), so the 
uncertainty of this estimate of 20 seconds is probably an order of magnitude. These 
soliton fronts are uncorrelated, although large scale structure should already be relevant 
even at z — 50. However it is hard to see how the birth pangs of different black holes could 
influence these fronts in a statistically detectable way. The round trip light travel time to 
the Moon is about 3 seconds, so a soliton passes through the Earth- Moon system roughly 
every ten round-trip travel times using the numerical estimates discussed above. The best 
way to detect the corresponding fluctuations would be to have an electromagnetic signal 
travel many times from the Earth to the Moon and back to see if there are transient shifts 
in an interference pattern every now and then (e.g. |91j). 



G. Ultra-precise timing 

In a very similar vein, ultra-precise timing experiments which test the steadiness of 
timing over time scales of order a few seconds to a few minutes should show these shifts 
as soon as the precision is high enough. Clearly the precision to detect such a signal has 
to correspond to seconds in the expansion rate of the universe, that is the precision must 
be of order 10~ 17 5 , or a few 10~ 18 . This precision is expected to be reached in the next 
generation of clocks [92H96] . 

Since a single soliton coming through has at least a transition time of the inverse of 
the frequency given above, the order of the transit time is some fraction of an entire 
day. New solitons come through on the order of magnitude of every 20 seconds, so how 
can we expect to be able to measure this? We note that a single soliton is a coherent 
superposition of approximately 10 97 correlated waves (using our adopted parameters, and 
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counting across the universe) combining to make up the soliton transition. The different 
solitons coming through every 20 seconds are incoherent. A careful error analysis might 
disentangle the two aspects. 

H. Summary of predictions 

This model makes quantitative predictions of several phenomena. Of the predictions 
discussed in this section, those in Sections C and G are the most likely to be verified or 
refuted in the near future. The predictions in the remaining sections are of cosmological 
importance, but they may require more time to test. 

VIII. CONSEQUENCES OF THE MODEL 
A. Primordial black holes 

In the approach described here primordial black holes grow - feeding on the background 

- as 

Mbh . / 3* 
mpi \t pi 

giving them a present mass of 3 x 10 16 g. 

In a speculative hypothesis in which all dark matter is black holes, a possibility which 
is not easily excluded (|97j). all primordial black holes would grow to 3 x 10 16 g. Therefore 
in our Galaxy the total number of primordial black holes is limited by 

Npbh 3 x 10 16 g < 1O 12 M ~ 10 45 - 3 g, (92) 

using the total mass of our Galaxy. If we hypothesize that all primordial black holes were 
initially at the Planck mass - an assumption from which it is straightforward to generalize 

- the relevant number of primordial black holes for our Galaxy is 

Npbh < 10 29 ' 3 , (93) 
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and their total mass is 



NpBHm Pl < 10 24 ' 3 g. 



(94) 



Comparing this to the total mass of the universe we see that 



N PBH mpi 
10 12 M m 



< 10 



-20.2 



(95) 



and therefore 



-20.7 



(96) 



If equality in this expression were true, then all primordial black holes would account for 
all dark matter. 

We note here an important hidden assumption: We implicitly assume that the primor- 
dial black holes do not get eaten by bigger black holes. Primordial black holes may be 
nearly at rest with respect to the microwave background system, and so may easily be 
eaten by the largest of the super-massive black holes. If that is an important process due 
to the evolution of large scale structure, then these limits derived here do not apply. 

There is a secondary difficulty in that under this hypothesis dark matter would not be 
constant, but would vary with cosmic epoch. Using the simple power-law approximation 
the mass dependence of the primordial black holes as a function of redshift would be 



which can probably be excluded already with existing data as the possible redshift de- 
pendence of dark matter. If this can be reliably proven, then the upper limit given above 
is firm. This limit is much stronger than discussed in |97j . 

We emphasize again, that black holes of low mass do not fulfill the occupation num- 
ber condition on spawning an active shell front. However, this condition just states 
according to the creation equation, that creation of new gravitons is strong, and per- 
haps self-sustaining. The creation equation clearly shows, as already emphasized above, 




(97) 
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that creation does continue, but at an extremely low level. Therefore, new gravitons are 
created all the time, but at a level energetically completely irrelevant compared to dark 
energy driven by the shell fronts from super-massive black holes. 

B. Lorentz Invariance 

It is clear from section IV.A that Lorentz Invariance is violated in our model in the 
limits of very high energy, a very large number of correlated bosons, very small length 
scales, and very small time scales. For instance, the minimum scale both in space and 
in time cannot be under-cut in any frame; similarly, the maximum energy cannot be 
exceeded in any frame. This is consistent with the Horava - Lifschitz theory of gravity 
in which Lorentz invariance is violated at very high energies or equivalently at very small 
distance scales [981 EH]. Even singularities cannot be pin-pointed with a precision better 
than these limits. This has led us to shells instead of horizons, and to a finite Lorentz 
factor for the shell front. We propose to deal with the further consequences elsewhere. 

C. Astronomical predictions 

In the scenario proposed here the main population of super-massive black holes has 
to be formed reasonably early to allow dark energy to be started. The efficiency of their 
contribution strongly decreases with redshift, and so dark energy is predicted to be very 
nearly constant to a fairly high redshift range, consistent with very recent work, |100j . 

Our approach strongly suggests that super-massive black holes grow mostly by merging, 
which should be testable with various observational strategies. 

Another consequence of our model clearly is, that reionization is extremely clumpy, 
and starts even before the first super-massive black holes, which in our model are formed 
out of the agglomeration of stars. Our model thus requires the early formation and 
demise of stars, and predicts the detection of signatures such as heavy elements, dust, 



53 



magnetic fields, and cosmic rays for such stellar events early in the history of the universe. 
Secondary consequences derive from cosmic ray interaction, such as high energy neutrinos; 
neutrinos should be strongest from the explosion of the first massive stars with sufficient 
heavy elements to form stellar winds to explode into (see, e.g., |101j ). 

D. Arrow of time 

In our model time is driven by the energy transfer from the Planck sea to our universe. 
Thus the arrow of time is uni-directional, since in our model energy flows in only one 
direction - from the background into our universe. An immediate consistency check of 
our model is the question of whether or not it describes a universe which is connected to 
the background Planck sea at all times. 

The horizon shells dip below criticality at the end of inflation, but at a much re- 
duced level they still create new gravitons, allowing a weak connection to the background 
throughout the transition phase from the cessation of inflation to the first creation of 
super-massive black holes. 

The shell fronts produced by the formation and mergers of the super-massive black 
holes come through at an estimated rate of one every twenty seconds, and any specific 
transition lasts a fraction of a day. So at any given moment we experience of order one 
thousand shell fronts, all connecting us to the background Planck sea. 

The shell fronts produced by stellar black holes or Planck particles do not give the 
full occupation number, and so their creation of new gravitons is enormously reduced, 
by many powers of ten. However, this creation is not zero. The number of stellar black 
holes is of order possibly 10 9 higher than the number of super-massive black holes, and 
their formation probably peaked relatively early in the universe, when the heavy element 
abundance was significantly lower (e.g. [90J). Therefore their repetition time is very much 
shorter (the time from one to the next), as is their transition time (the step-length, which 
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is essentially the inverse gravitational wave frequency). So over a very wide range of 
time-scales there is always a connection to the background, albeit extremely weak. So the 
arrow of time is always assured. 

In the far distant future, when the shell fronts of the super-massive black holes be- 
come sub-critical, causality will always be assured by the decreasing interaction with the 
background. The universe will no longer be driven by dark energy, as it will turn into a rel- 
ativistic gas, expanding essentially adiabatically. The universe will end altogether, when 
this tenuous connection to the background lapses, since causality is no longer possible. 



E. Signal speeds and mass 

The Lorentz factor for the effective propagation of a signal derived above should apply 
to the propagation of any signal carried by mass-less boson particles, or waves travelling 
with the speed of light even at distances as small as sub-atomic scales. In our model such 
signals never travel at exactly the speed of light, but always just very slightly below. The 



limit is given by the number of coherent waves (see Eq. 49 ) 




where E\ is the energy of a single boson, E^, is the total energy of all coherent bosons, 
and d is the distance travelled. The smallest possible limit is obviously the Planck length, 
and at that limit the Lorentz factor becomes 

d 



r 



(99) 



4/«. 

At the smallest distances sometimes considered in particle physics, of order ~ lCr 18 cm, 
this Lorentz factor would still be T r ~ 10 +7 ' 5 , but would require, independent of the 
length d, that 

E 1 E^ = 2(2vr) 3 (m Pl c 2 ) 2 ~ 10 35 - 3 erg 2 . (100) 
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It is straight-forward to verify that shell fronts obey this condition. 

However, there is a limit to this kind of reasoning, since the smallest possible energy 
unit in the universe is the energy of a gravitational wave spanning the entire universe. 
This energy is 

ch 

Ex = ~ l(T 43 - 9 erg. (101) 

If the radius of the universe is taken to be a half wavelength, then the length scale is twice 
as large and the energy is half as large. This energy is the equivalent energy obtained by 
using the Gibbons-Hawking temperature ([21 E])- At high redshift this scales as (1 + z) 3 / 2 , 
but with a factor 1.3 due to the subtle effect of dark energy. At the end of inflation a 
phase transition occurs. Requiring that the phase transition parameter e = 1 in natural 
units, the energy E\ is 10~ 8 erg which implies that (z in fi ~ 10 29 ) 

E^tnfl- 100GeV(^) 3/2 , (102) 

as the energy where mass first begins to manifest itself in the evolution of the universe. 
When the first rest mass appears, the scaling of energies with look-back time fails and 
inflation stops. This value is obviously subject to the uncertainty of when inflation really 



ends. This can also be written (see Eq. 84) to within a factor of order unity as 

E ltinfl ~ m Pl c 2 e -W^ . (103) 
Integrating over all gravitons created during inflation by number and energy shows that 



Eq. |100| is obeyed right through the inflationary period, with E\ always the lowest energy 
in the system. 

This is a crude estimate of the individual mass of a boson at the end of inflation, 
where we are able to reach the Planck uncertainty, imposing a behavior that first gives 
the outside appearance of mass. This energy is of the same order of magnitude as that of 
the Higgs boson (see FigJT]). 
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FIG. 7: The lowest energy graviton existing today is the redshifted lowest energy graviton at 
the end of inflation. At the end of inflation the energy of this graviton was ~ 100 GeV. 

In a future investigation we will carry out a more refined mathematical calculation to 
determine whether or not there is a connection between the signal mass at the end of 
inflation and that of the Higgs boson. 

This is the signature of a finite rest-mass. The fact that wave pulses travel at very 
slightly subluminal velocities is due, in the large number limit, to the uncertainties in 
space and time at the Planck level. 

F. Consistency 

The model of the generation of gravitons by shells (shell fronts, Planck shells, or horizon 
shells) is not a fully developed theory. What we have shown is that an approach exists, 
which given its assumptions can provide a common basis for the understanding of such 
diverse phenomena as black holes, dark energy, inflation, rest mass, and the arrow of time. 
Clearly this is merely a first step in the development of a complete theory. 
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A fully developed theory would start with some prescription for the production of 
various black holes and their interactions including mergers. A complete theory could 
describe their formation, growth and merger history (e.g. [351 G2E UU H02H111] ). 
and at the same time describe the evolution of their shell fronts, which contribute to dark 
energy. Over a certain mass range black holes drop out of creating dark energy, since they 
can no longer provide the critical surface density on their shell front. The description of 
all these processes can be given in terms of functions of the redshift z. Dark energy 
enters the integrals for luminosity distance and Hubble parameter evolution, turning the 
creation equation into a set of coupled non-linear differential and integral equations. The 
overall evolution of the universe is then governed by the creation equation integrated over 
all of the black holes, over their creation and over their merger history. 

In the three applications of the creation equation we used gravitons pulled from the 
background usually below the maximum energy in the background, so that with the proper 
Lorentz factor they are at the Planck peak on the shell. The actual number pulled is far 
below the density in the background for the two applications of the shell front and the 
black holes. But only for the horizon shell, so for the inflationary phase, the density of 
gravitons pulled in is so high, that modifications of the background spectrum due to the 
appearance of separate forces of nature do become important. Only in that case do we 
reach the limits of the phase space density in the background itself. 

The model we propose in this paper could be disproved in several ways: a) if the 
gravitational wave background corresponding to dark energy is not detected with the 
appropriate strength and spectrum; b) if the formation of the super-massive black holes 
is near redshift 10 or below rather than near 30 - 50; or c) if the ultra-precise timing and 
lunar laser ranging does not detect the individual soliton fronts of gravitational waves. 
If, however, these predictions were confirmed, new paths of investigation would open up 
for research on dark energy, black holes, inflation, Lorentz transformations, the arrow of 
time, and the characteristic energy at which rest mass appears. 
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G. Summary 

We have introduced the notion of stimulated emission of gravitons on a critical surface, 
drawing energy from a background field. The phase space density used in our model is 
derived in close analogy to the holographic principle. The equation derived in Section 
III for the stimulated emission of gravitons from a background is the first new element 
in our approach, and the use the equation for the very large Lorentz factor derived from 
Planck uncertainties is the second new element. This allows the interpretation of dark 
energy, black holes, the inflationary period of the universe, the arrow of time and the 
characteristic energy at which rest mass first manifests itself to be put on a common 
basis. The other concepts used in our model are variants of commonly discussed ideas in 
gravity theory. 

The next step is to develop a more formal mathematical theory to self-consistently 
derive the creation equation. Part of this formal development will include a proper treat- 
ment of Lorentz transformations. An additional step is to include the spin of the black 
holes, and another is to explore some of the consequences of our model in more detail, 
such as the signal speed limit; the signal speed limit directly relates to the characteristic 
energy at which rest mass first appears. Most importantly, we will develop more specific 
proposals for experimental tests of this general approach. 
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IX. APPENDICES 

A. Graviton distribution function 

The graviton distribution function on the weak-gravity brane (our brane) for a Planck 
surface interacting with a background field can be determined from the Boltzmann equa- 
tion for coherent production of gravitons in a Friedman- Robertson- Walker (FRW) universe 
[112] . see [113] . The distribution function J\f(k,t) satisfies the equation (h = c = 1 for 
the derivation of the differential equation satisfied by Af(k,t)) 

( d _k(t) d\ i r d 3 k' r d 3 P r d 3 P ' r d 3 k" 

[di ~ Btfj dk) ^ '*) - l J (2tt)3 2A;' J {2nf2E{p) J (2vr) 3 2 E{p>) J (2tt) 3 2 k" 

\M\ 2 2 7c 4 5\k +p-k' -k" - p') 2 vr 3 5 3 {k> - k") 
(M(k f , t) M(k", t) g b (p', t) (1 + J\f(k, t)) 
- Af(k, t) g b (p, t) (1 + Af(k', t)) (1 + Af(k", t))) (104) 

The 5-function, 5 3 (k' — k"), has been inserted to impose coherence of the outgoing gravi- 
tons. |M| 2 is the matrix element squared for the quadrupole emission of a graviton of 
energy k". gb{p,t) is the occupation number distribution of the background particle sea. 
R(t) is the scale factor for an expanding universe. 
The 4-dimensional 5-function can written as 

5\k +p-k'-k"- p') = 5{k + E p -k' - kl - E' p ) 5 3 (k +p-k'- k" - p>) (105) 
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Imposing the coherence condition and assuming that the momentum transfer k — k! is 
small, the <5 3 -term can be expanded in a series 

5 3 (k + p~2k' -p') = 5 3 (p-p'-k') + (k-k')—5 3 (p-p'-k') 

+ \(Ak- ^5 3 (P- P' -#) + ••• (106) 
where q = p' + k' . If the background particles are assumed to have structure, to have 

— * 

excited states and to obey the relation p ~ p', then first term in the expansion of 
5 4 (k + p — k' — k" — p') can be factored as 

5\k +p-k' -k" - p') = 5 3 (p- p' - k>) 5{k - k' ) (107) 

— * — * 

where the coherence condition k' = k" has been inserted. The second term in the expan- 
sion vanishes when it is integrated over all angles. 

The Boltzmann equation for Af(k, t) to lowest order in the expansion of the 4- 
dimensional 5-function is 

B-W) k i) mj ^^ mj)(mj)+l) ' (108) 

where the factor n is given by 

J_ f d 3 p [Ml 2 
2nR(t)J (2tt) 3 j B(p) 2 " R(t) 1 ' 

and |M| 2 , as stated above, is the matrix element squared for quadrupole emission of a 

graviton wave from the background. 

The product Af(k,t) (J\f(k : t) + 1) can be expressed as 

- Af(k, t) (Af(k, t) + l) = ^T a ^Af(k, t) , (110) 

where T g = T g0 /R(t) and k is the magnitude of the graviton momentum. After making 
the replacement k = k/R(t), the equation for M(k,t) becomes 

!^M = «o^i^(M), an) 

ot c kok 
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or in terms of the frequency of the wave at emission vq 

The solution of this equation can be written in terms of dimensionless parameters x and 
y, where 

Kb T 



and 

rt K ° C T9 °dt' (114) 



k B T T 



as 

i x 3 
ir' 2 e x — 1 

with 



tf(x,y) = —- 5r ^ (115) 



X = y/{x* + 2y} . (116) 

When y reaches the same order of magnitude as x 2 , the number density Af(x, y) begins to 
deviate from the Bose-Einstein distribution. This corresponds to the appearance of the 
Higgs boson. 



B. The Kompaneets equation 

For photons, also bosons, the basic equation is the Boltzmann equation describing the 
scattering of photons on electrons, and includes the differential cross-section, an integral 
over directions, the distribution function of the electrons, and the phase space distribution 
of the photons, always including stimulated emission which implies the use of the square 
of phase space density, see above, and also, e.g. [S]. In this application the equation was 
first derived by Kompaneets. 
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The relativistic Boltzmann equation describing stimulated emission has already been 
given above. The non-relativistic form for stimulated emission for photons scattering from 
a distribution of electrons 



-n{u) = c d'p —dfi(/ e (p 1 )n(a; 1 ){l + n(a;)}-/ e (p)nH{l+n(a; 1 )})(118) 



where f e {p) is the distribution of electrons, n(u) is the distribution function of photons, a 
is the cross-section, and the photon density squared term accounts for stimulated emission. 

Simplifying this Boltzmann equation to the limit that the energy change is minimal, 
that the cross-section is isotropic, the phase space density n(u) is very large, and the 
distribution functions are very peaked, we obtain 



neglecting lower terms, which do not scale with n 2 (u)). Here A{pc} is the energy being 
exchanged. We will use this expression for insight as to how an equation may be structured 
to describe stimulated emission of gravitons from a background. Of course this is just 
a didactic step, since the production of gravitons in our world by stimulated emission 
from the background is clearly relativistic. What we want to describe is the total energy 
injection by stimulated emission, so we have to modify this expression to account for the 
energy emitted into our world, and also to account for the total phase space affected in a 
given time interval. 



p + u <=> p 1 + Ui 



(117) 



is (see, e.g. [9]) 





(119) 



C. 



The end of dark energy 



Dark energy runs out of strong creation when the background phase space density 
becomes too low at resonance to supply enough gravitons to keep dark energy production 
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going. The parameter y as defined in Eq. 114 reaches the same order of magnitude as x 2 
when 



x 2 



y ~ — . (120) 
The future scale factor of the universe is 



R(t) = R e Hot (121) 



The parameter y reaches the scale of x 2 at time t^E, which scales as 



t DE = ^{-ln(l + z)}. (122) 



and is determined from the relation 



x 2 ~ 10-» 2 (l + z DE ) 2 ~ ^gt DE (123) 



or in Planck units 



J^e^-p* (121) 



Solving for z^jg implies that e, the dimensionless value of Ko, is of order (tpiH ) 2 . A given 
value for e determines the time at which dark energy will cease to be renewed. Thereafter 
the universe will expand adiabatically. 



D. The cosmological equations 

The effect of dark energy on the metric tensor elements which describe the evolution 
of the universe can be determined by calculating the energy-momentum tensor for a fluid 
whose entropy continually increases. The metric for our universe is taken to be of the 
standard form as shown in the invariant line element 



ds 2 = -c 2 dt 2 + a(t) 2 (dx 2 + dy 2 + dz 2 ) 



(125) 
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A Lagrangian density of the form [114] 

1 r 



VT 



-9 (p + Pde) (c 2 + S) 



(126) 



where £ is the rest, internal specific energy of the fluid, leads to the desired form of the 
energy- momentum tensor, provided that, 



d£ _P 

dp p 2 



(127) 



and 



d£ 



DE 



dpDE P 2 

The total energy momentum tensor is assumed to be the sum of two contributions 



128) 



rpfiV rp(F)fJ,U _|_ rp(DE) 



(129) 



where T^ u and T^ DE ^ U are the energy- momentum tensors for a fluid in our universe and 
for the dark energy in our world, extracted from the background Planck sea, respectively. 
Assuming that the dark energy components, T^ DE ^ U , are functions of t only, T^ u can be 
written as 



rpflU 



p(t,x,y,z)+T°° E (t) 



U {a(t)f + 1 DE l r J 














P(x,y,z) 

Ht)) 2 






+ n% w 









pi P{x,y,z) , rp33 ,.\ 

U (a(t)) 2 DE W 



(130) 



The conservation of energy relations are 



T"".^ = 



(131) 
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rp/j.0 (f tP (t,x,y,z ))a(t) + (iTOQ E (t))a(t) _ 3 (|a(i)) p (x,y, z) + (|a(t)) T& (t) (a(t)) 2 



+ 



a(t) 



+ 



a(t) 



3 (|a (Q) p (t, x, y, g) + 3 (|a (t)) rg| (f) + (|a (t)) 

a(t) 

, r^(t)(a(t)) 2 + (^o(t))r« (t)(o(t)) 2 

a(t) 



(132) 



(a(t)r 
(a(i)) 2 



(a(t)Y 



The last three equations can be written as 



(133) 



VP = 

If the dark energy is considered to be a kind of fluid, then T^ v has the form 

p(t,x,y,z) + g{t)p DE (t) 















(134) 










(135) 



and the covariant divergence equations become 

(|p (t, x, y, z)) a (t) + {f t g(t)p DE (*)) a (t) + 3 P (x, y, z) |a 



rpfj.0 

if* 



a(t) 

+3 (£a{t))f{t)P DE (t) (a(t)) 2 
a(f) 

+3 (|a (Q) p (t, x, y, *) + 3 (|a (t)) (t) 

a(i) 



(136) 
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Conservation of energy requires = 0. In order for the numerical factor and the 

powers of the scale parameter a(t) to match we have to identify the central terms 

d v , 3 (ia(t))f(t)P DE (t)(a(t)) 2 3(ia(t))g(t)p DE (t) 

it mpDE) + W) W) = ( } 



or 



(g(t)p DE ) + 3 Q (f(t) a 2 P DE + g(t) p DE ) = (138) 
The source term in the conservation equation is 3po E H(z). This term can also be 



d 
dt 



obtained by first combining Eqs.50 and 53 with the density, ubh,o, of black holes to 
obtain 

Pde = ^ M BH c 2 riBHfl (139) 



and then using [54] to obtain the rate of dark energy production 

£de = g M bh c 2 n BH ,o H{z) = 3 p DE H(z) (140) 

The expressions for g(t) and f(t) are not fixed by either Einstein's equations nor the 
continuity equation. Therefore we set g(t) = 1 and determine a reasonable form for 
f(t) by the following argument. Gravitational systems may have a negative specific heat 
[115] , and examples are Schwarzschild black holes, and self-gravitating globular clusters of 
stars. If BHs are in a non-flat geometry such as de Sitter or anti-de Sitter (cosmological 
constant different from zero), then there are values for the angular momentum (Kerr 
BH) and charge (Reissner-Nordstrom BH) for which the specific heat is negative, but 
these are just examples. This is exactly the situation we are proposing. So we suggest 
that f(t) = —1/a 2 , or equivalently that Pd E = —Pdec 2 . This solves, as usual, the 
requirement that normal matter runs as (1 + z) 3 , and that the conservation of energy is 
obeyed. 
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The Einstein equations are automatically satisfied, since we use the normal energy- 
momentum tensor. With these choices for g(t) and f(t) the equation satisfied by a(t) can 
now be determined. The Einstein equations can be written as 

iV = -8nG N (t^ - l - 9iiv Tj (141) 

where T = T x A . The t, t component of the Ricci tensor is 

t^a (t) , . 

Ru = S^-^ 142 
a {t) 

and the spatial components of the Ricci tensor all have the form 

/ d \ 2 d 2 

Rii = ~ 2 (jt a ^) ~ a ®di? a ® 2 = 1 > 2 ' 3 - ( 143 ) 

These two equations and the expressions for the can be used to eliminate the a(t) 
term, giving 

f da(t)\ 2 87rGjva(i) 2 



V dt ) 



(p(t)+g(t)p DE (t)) (144) 



3 

This is expression for (da/dt) 2 is of the standard form, but with p(t) replaced by 
pit) + g(t) pDE{t)'i as shown above the most likely solution is that g(t) = 1. The only 
modification which we have made is to identify the dark energy instead of subsuming it 
into the density as one term. 

This completes this argument; our proposed theory for the origin of dark energy satisfies 
the energy conservation equation including the proper source term, and the Einstein 
equations. 



E. Graviton emission from black holes 



Here we derive the steady graviton emission from black holes in a manner analogous 
to Parikh & Wilczek Q55J) in an expanding universe. 
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Defining 

p = r e m/2 , and m = /x(0) = fi(t) e m/2 ; e m/2 = a(t) , (145) 



Eq. ( 70 ) can be written as 

ds 2 = ( ; + ^t 2 + (l + ^) [V + P 2 (146) 

The emission of a graviton of energy E g = hu lowers the mass of the black hole by the 
amount E g to m — u (h = c = 1). So the metric after emission is 

The action for an outgoing S-wave particle is 

"Pout 



rpout 

S = / ppdp 

J Pin 



pout rm{t)—w(t) 

- dp (148) 



' Pin Jm(t) P 

where if = m(t) — uj(t), p in = a(t)m, p out = m{t)—u{t). The null radial p geodesies 
yield 



dp _ 1 ~ 2 P 

dt ( , (-m(t)-u;(f)) N 3 



(149) 



1 ' 

Inserting this expression for the p term, the action becomes 

rPout r^fl + ^mY 

s = - / \ (m(t) LJ dJd p ( 15 °) 

•J Pin JO 1 Y~ p 

Performing the p integration first and then using contour integration to carry out the 
integration over u', the imaginary part of the action for an S-wave particle crossing the 
horizon in an outward direction is found to be 

= 2na(t) 2 uj(hm- |) (151) 
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F. List of definitions 

1. "Planck" and "shell" language 

We use many words with the prefix "Planck" , most of them are well established, and 
only very few of which are new; we also use "shell" several times: 
Planck length: l 2 Pl = G N h/c 3 , with l Pl ~ 1.6 • l(T 33 cm. 
Planck time: r 2 Pl = G N h/c 5 , with r Pl ~ 5.4 • 10~ 44 s. 
Planck area: a P i = 7tGat/i/c 3 , with a P i ~ 8.2 • lCT 66 cm 2 . 
Planck mass: m Pl = hc/G N , with m P i ~ 2.2 • 10~ 5 g. 
Planck energy: E Pl = hc b /G N , with m P{ ~ 2.0 • 10 16 erg. 
Planck spectrum: The standard Planck spectrum for even spin particles. 
Planck sea: The background field, assumed here to have a Planck spectrum at maximal 

density and energy. 
Planck shell : The surface corresponding to a black hole. 

Shell front: The surface of the graviton shell flying out from the formation of a black 

hole, or the merger of two black holes. 
Horizon shell : The surface of interaction with the background at the horizon during 

inflation. 

2. Redshift Formulation of Time and Distance 

In order to evaluate the behavior of gravitons we need to use integrals which allow 
an extrapolation into the future, implying (1 + z) — > 0. We set the low redshift to z, 
and the high redshift to z±, the redshift at which we assume the main gravitons to be 
produced. 



70 



The expression for the time can be written the normal way 



dz' 



(152) 



t = 



(1 + z>) H(z') ' 



where 1 + z can have any value above zero. The time into the future diverges with the 
natural logarithm of the redshift expression 1 + z going to zero. The derivative of redshift 
z with respect to time can be written as 



The corresponding integral for the luminosity distance between redshift z and redshift 
z* can be written as 



and again this is valid into the future. We use redshift z± as the redshift, when a black 
hole was formed, or a merger between two black holes happened, so in our interpretation 
the epoch, when the shell front was initiated in the last few Planck times before what is 
normally called horizon formation. Similarly the luminosity distance diverges as (1 + z)^ 1 
in the future. 



[1] Grishchuk, L. P., Class. & Quant. Grav. 10, 2449 (1993). 

[2] Gibbons, G. W., Hawking, S. W., Phys. Rev. D 15, 2752 (1977). 

[3] Gibbons, G. W., & Hawking, S. W., Phys. Rev. D 15, 2738 (1977). 

[4] Harms, B., & Leblanc, Y., Phys. Rev. D 46, 2334 (1992). 

[5] Harms, B., & Leblanc, Y., Phys. Rev. D 47, 2438 (1993). 

[6] Bekenstein, J. D., Phys. Lett. B 481, 339 (2000). 

[7] Bousso, R., Rev. Mod. Phys. 74, 825 (2002). 

[8] Kompaneets, A.S., Sov. Phys. J. Exp. & Theor. Phys. 4, 730 (1957). 



dz 
dt 



(l + z)H(z). 



(153) 




(154) 



71 



[9] Rybicki, G.B., & Lightman, A. P., Radiation processes in Astrophysics, Wiley (1985). 

[10] Zel'dovich, Ya. B., Usp. Fix. Nauk 115, 161 (1975); transl. Sov. Phys. - Usp. 18, 79 (1975). 

[11] Frieman, J. A., Turner, M. S., & Huterer, D., Annual Rev. of Astron. & Astrophys. 46, 

385 (2008). 

[12] Lemaitre, G., Nature 127, 706 (1931). 

[13] Binetruy, P., et al. Phys. Lett. B 477, 285 (2000). 

[14] Randall, L., Sundrum, R., Phys. Rev. Letters 83, 3370 (1999a) 

[15] Randall, L., Sundrum, R., Phys. Rev. Letters 83, 4690 (1999b) 

[16] Dirac, P., Nature 139, 323 (1937). 

[17] Zel'dovich, Ya. B., JETP Letts. 6, 316 (1967). 

[18] Weinberg, S., Rev. Mod. Phys. 61, 1 (1989). 

[19] Bekenstein, J. D., Phys. Rev. D 9, 3292 (1974). 

[20] Bekenstein, J. D., Phys. Rev. D 12, 3077 (1975). 

[21] Bekenstein, J. D., Phys. Rev. D 23, 287 (1981). 

[22] Bekenstein, J. D., Phys. Rev. D 27, 2262 (1983). 

[23] Bekenstein, J. D., Phys. Rev. D 49, 1912 (1994). 

[24] Bekenstein, J. D., Contemp. Phys. 45, 31 (2004). 

[25] Macdonald, D., & Thorne, K.S., Month. Not. Roy. Astr. Soc. 198, 345 (1982). 

[26] Macdonald, D.A., & Suen, W.-M., Phys. Rev. D 32, 848 (1985). 

[27] Price, R.H., & Thorne, K.S., Phys. Rev. D 33, 915 (1986). 

[28] Thorne, K.S., k Macdonald, D., 1982 Month. Not. Roy. Astr. Soc. 198, 339 (1982). 

[29] Feynman, R.P., Leighton, R.B., Sands, M., The Feynman lectures on physics: Quantum 

mechanics, 1965, Addison- Wesley, Reading, MA, USA 

[30] Higgs, P.W., Phys. Rev. Letters 13, 5081 (1964); Phys. Rev. 145, 1156 (1964). 

[31] Englert, F. and Brout, R., Phys. Rev. Letters 13, 321 (1964). 

[32] Hawking, S. W., Phys. Rev. Letters 26, 1344 (1971). 



72 



[33] Wald, R. M., Phys. Rev. D 56, 6467 (1997). 

[34] Caramete, L.I., & Biermann, P.L., Astron. & Astroph. 521, id.A55 (2010), 

larXiv:0908.2764l 
[35] Berti, E., Volonteri, M., Astrophys. J. 684, 822 (2008). 
[36] Gergely L. A., Biermann, P.L., Astrophys. J. 697, 1621 (2009). 



[37] Biermann, P.L., Kusenko, A., Phys. Rev. Letters 96, 091301 (2006); |astn>ph /0601004. 
[38] Silk, J., Tracing the Ancestry of Galaxies (on the land of our ancestors), Proa, IAU 

Symposium, 27, 273 (2011). 
[39] Hawking, S. W., Phys. Rev. D 13, 191 (1976). 
[40] Bekenstein, J. D., Phys. Rev. D 7, 2333 (1973). 

[41] Bardeen, J. M., Carter, B., Hawking, S. W., Comm. Math. Phys. 31, 161 (1973). 
[42] Cook, G. B., & Abrahams, A. M., Phys. Rev. D 46, 702 (1992). 

[43] Joshi, P. S., & Malafarina, D., Int.J.Mod.Phys. D20, 2641(2011); eprint ' arXiv: 1201 .36601 
(2012). 

[44] Courant, R., Friedrichs, K.O., Supersonic flow and shock waves, Springer- Verlag (1948), 

reprinted in its entirety (1976). 
[45] Susskind, L., & Lindesay, J., An introduction to black holes, information, and the string 

theory revolution : the holographic universe, World Scientific, (2005). 



[46] Almheiri, A., et al. eprint arXiv: 1207.3123 (2012) 



[47] 't Hooft, G., Phys. Scr. T 36, 247 (1991), Proc. conf. The birth and early evolution of our 

universe; Proceedings of Nobel Symposium 79. 
[48] Hawking, S. W., Comm. Math. Phys. 43, 199 (1975). 
[49] Page, D. N., Phys. Rev. D 13, 198 (1976a). 
[50] Page, D. N., Phys. Rev. D 14, 3260 (1976b). 
[51] Page, D. N., Phys. Rev. D 16, 2402 (1977). 
[52] Hawking, S., Month. Not. Roy. Astr. Soc. 152, 75 (1971). 



73 



[53] Carr, B. J., Astrophys. J. 201, 1 (1975). 

[54] Sanchez, N., Phys. Rev. 18, 1030 (1978). 

[55] Parikh, M. K., Wilczek, F., Phys. Rev. Letters 85, 5042 (2000). 

[56] McVittie, G.C., Month. Not. Roy. Astr. Soc. 93, 325 (1933). 

[57] LIGO and Virgo collaborations, Nature 460, 990 (2009). 

[58] Shakura, N. I., & Sunyaev, R. A., Astron. & Astroph. 24, 337 (1973). 

[59] Novikov, I. D., & Thorne, K. S., in Black holes (Les astres occlus), Eds. DeWitt C, DeWitt 

B., Gordon & Breach, New York, p. 343 (1973). 

[60] Bardeen, J.M., Nature 226, 64 (1970). 

[61] Guth, A. H., Phys. Rev. D 23, 347 (1981). 

[62] Linde, A. D., Phys. Lett. B 99, 391 (1981). 

[63] Linde, A. D., Phys. Lett. B 108, 389 (1982). 

[64] Sanchez, N., Latest News from the Universe, Torino Cosmology Colloquium 
held 21-24 Oct., 2009 at Torino Science Academy of Turin, Torino, Italy. At 



http : / / ecolechalongetorino.oato. inaf . it , p . 1 . 



[65] de Vega, H. J., & Sanchez, N. G., eprint |ailGv:astro-ph/070m^ (2007). 



[66] Stanev, T., High energy cosmic rays, Springer-Praxis books in astrophysics and astronomy. 
Springer, Heidelberg (2010), new edition. 



[67] Greene, J.E., Barth, A.J., k Ho, L.C., New Astron. Rev. 50, 739 (2006); arXiv:astro- 
|ph/05li810] 

Greene, J.E., & Ho, L.C., Astrophys. J. 670, 92 (2007a). 
[69] Greene, J. E., & Ho, L. C, in proc. The Central Engine of Active Galactic Nuclei, ASP 
Conf. 373, 33 (2007b). 



[70] Greene, J.E., Ho, L.C., & Barth, A.J., Astrophys. J. 688, 159 (2008); |arXr^0 810. 1972 
[71] Gilmore, G., et al., Astrophys. J. 663, 948 (2007). 

[72] de Vega, H. J., & Sanchez, N. G., Month. Not. Roy. Astr. Soc. 404, 885 (2010). 



74 



[73] Portegies Zwart, S. F., et al., Nature 428, 724 (2004). 

[74] Sanders, R. H., Astrophys. J. 162, 791 (1970). 

[75] Cote, P., Piatek, S., Ferrarese, L. et al. Astrophys. J. Suppl. 165, 57 (2006). 

[76] Balcells, M., et al. Astrophys. J. 665, 1084 (2007). 

[77] Graham, A. W., Month. Not. Roy. Astr. Soc. 422, 1586 (2012). 

[78] Scott, N., Graham, A. W., eprint |arXiv:1205.5338| (2012). 



[79] Leigh, Nathan, Boker, Torsten, Knigge, Christian, Month. Not. Roy. Astr. Soc. 424, 2130 



(2012); arXiv: 1205.5033, 



[80] Neumayer, N., Walcher, C. J., Adv. in Astron. 2012, id. 709038 (2012). 
[81] Yungelson, L.R., et al., Astron. & Astroph. 477, 223 (2008). 
[82] Appenzeller, I., & Fricke, K., Astron. & Astroph. 18, 10 (1972a). 
[83] Appenzeller, I., & Fricke, K., Astron. & Astroph. 21, 285 (1972b). 
[84] Allen, B., Phys. Rev. D 37, 2078 (1988). 

[85] Sesana, A., Vecchio, A., k Colacino, C. N., Month. Not. Roy. Astr. Soc. 390, 192 (2008). 
[86] Sesana, A., & Vecchio, A., Phys. Rev. D 81, id. 104008 (2010). 

[87] Kramer, M., in Relativity in Fundamental Astronomy: Dynamics, Reference Frames, and 
Data Analysis, Proc. IAU Symp. 261, p. 366 (2010). 
Liu, K. et al, Month. Not. Roy. Astr. Soc. 417, 2916 (2011). 



[89] Wyithe, J.S.B., Loeb, A., Astrophys. J. 590, 691 (2003); arXiv:astro-ph/0211556 



[90] Mirabel, I.F., et al., Astron. & Astroph. 528, id.A149 (2011). 



[91] Steinmetz, T., et al., Science 321, 1335 (2008); arXiv:0 809.1663 

[92] Chou, C. W., et al., Science 329, 1630 (2010). 

[93] Chou, C.W., et al., Phys. Rev. Letters 104, 070802 (2010). 

[94] Parthey, Ch. G., et al., Phys. Rev. Letters 107, 203001 (2011). 

[95] Predehl, K., et al., Science336, 441 (2012). 

[96] Diddams, S. A., et al., Science 306, 1318 (2004). 



75 



[97] Abramowicz, M. A., et al., Astrophys. J. 705, 659 (2009). 

[98] Hofava,P., Phys. Rev. D D 79, 084008 (2009). 

[99] Visser, M., J.Phys.Conf.Ser. 314, 012002 (2011); arXiv:1103:5587. 



[100] Reid, B. A., et al., Month. Not. Roy. Astr. Soc. , 426, 2719 (2012); eprint arXiv:1203.6641 



[101] Biermann, P.L., & Cassinelli, J.P., Astron. & Astroph. 277, 691 (1993); astro-ph/9305003 

[102] Abel, T., et al., Science 295, 93 (2002). 

[103] Bellovary, J., et al., Astrophys. J. 742, id.13 (2011). 

[104] Fanidakis, N., et al, Month. Not. Roy. Astr. Soc. 419, 2797 (2012). 

[105] Li, Y.-R., et al., Astrophys. J. 742, 33 (2011). 

[106] McMillan, S. L. W., Portegies Zwart, S. F., in Proc. Massive Stars in Interactive Binaries, 
ASP Conf. Ser. 367. Edited by N. St.-Louis & A. F.J. Moffat. San Francisco: A.S.P., p.697 
(2007). 

[107] Portegies Zwart, S. F., & van den Heuvel, E. P. J., Nature 450, 388 (2007). 

[108] Portegies Zwart, S. F., McMillan, St. L. W., & Gieles, M., Annual Rev. of Astron. & 

Astrophys. 48, 431 (2010). 
[109] Volonteri, M., et al., Astrophys. J. 582, 559 (2003). 
[110] Volonteri, M., et al., Month. Not. Roy. Astr. Soc. 383, 1079 (2008). 
[Ill] Volonteri, M., et al., Month. Not. Roy. Astr. Soc. 384, 1387 (2008). 

[112] Weinberg, S., Gravitation and Cosmology: Principles and Applications of the General 

Theory of Relativity, Wiley- VCH, 688 pp., (1972). 
[113] Bernstein, J., & Dodelson, S., Phys. Rev. D 41, 354 (1990). 
[114] Ray, J., J. Math. Phys. 13, 1451 (1972). 
[115] Spitzer, L., Jr., Astrophys. J. Letters 158, L139 (1969). 



